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ABSTRACT 


The  aim  of  this  research  was  to  obtain  a quantitative  understanding  of  the 
McHugh  stean  water  pulsejet  cycle,  which  in  its  simplest  embodiment,  is  a 
thrust-producing  engine  with  no  moving  parts.  The  cycle  is  also  adaptable 
to  water  pumping,  a study  of  which  is  currently  being  funded  by  I'.RDA;  or, 
indeed,  to  the  pumping  of  any  vaporizable  fluid.  It  has  also  i)een  used  as 
an  agitator,  and  may  find  applications  as  such  where  a fluid  or  slurry  is 
too  corrosive  for  conventional  mechanically-driven  pumps  or  agitators.  The 
basic  theory  may  also  be  peripherally  helpful  in  such  diverse  fields  as  the 
"chugging"  of  atomic  reactors  during  emergency  shut-down  and  the  catastrophic 
explosions  which  can  result  from  the  dynamic  mixing  of  v^ater  with  lava  or 
molten  metal. 

Although  by  no  means  complete,  it’s  believed  that  the  theory  presented  herein 
adequately  explains  the  McHugh  cycle,  and  points  the  way  for  further  perfor- 
mance improvements.  Steam  water  pulsejets  arc  not  yet  as  efficient  as  con- 
ventional steam  engines,  but  tlierc  may  be  applications  V'jhcre  the  rather 
extreme  mechanical  simplicity  makes  them  cost  effective.  In  terms  of  "speci- 
fic fuel  consumption,"  the  best  engine  tested  cori'csponded  to  about  0.56  lb 
of  fuel  per  hour  per  pound  of  thrust,  assuming  an  80^  boiler  efficiency.  This 
is  comparable  to  a turbojet,  but,  of  course,  the  pulsejet  has  a natural  ad-  _ i 
vantage  in  the  <Jenscr  medium.  TMrs'  fuel  consumption  is  about  five  times 
we- would  expect  'of  a diesel  engine  driving  a water  propeller,  p^ut  it  is 
comparable  with  what  one  might  expect  from  a fractional  horsepower  I.C.  engine 
driving  a water  propeller;  which  is,  after  all,  the  direct  comparison  in  scale 
for  the  laboratory  engines  so  far  tested.  Since  the  technology  is  new, 
one  would  expect  improved  results  with  further  work. 

The  pulsejets  tested  discliarged  essentially  all  the  boiler  heat  to  the  test 
tank.  The  longest  one  tested  was  7.75  feet,  but  there  was  no  indication  that 
this  was  an  upper  limit.  (Performance  generally  improves  somewhat  with 
increasing  length.)  The  highest  heating  rate  achieved  was  10  Kw,  and  for  a 
0.875  inch  I.D.  duct,  this  represents  a larger  heat  flow  than  is  achievable 
with  a conventional  water  filled  licat  pipe.  Thus,  the  technology  may  find 
applications,  in  addition  to  those  already  mentioned,  in  the  transfer  of  heat. 
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This  stuJ>’  tsT'.'jraco.s  disciplines  which  nrc  rarely  brought  together  in  a single 
analysis.  The  writer  therefore  had  the  choice  of  using  entirely  unfamiliar 
symbols  for  nany  variables,  in  order  to  avoid  duplication,  or  using  the  usual 
ones.  The  latter  course  was  cltosen,  so  that  (for  exami)lcj  at  various  places 
Cp  is  a power  coefficient,  a specific  heat  or  a pressure  coefficient.  It's 
believed  that  no  confusion  can  arise  from  this;  where  confusion  could 
arise,  some  symbols  have  been  changed. 
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INTRODUCTION 


This  I'.ipcr  is  concerned  with  a remarkable  fluid  flow  phenomenon  which  should 
lo};icail/  ho  called  ''the  Mclliiyh  effect,"  oi'  "McHugh  cycle"  after  its  discov- 
erer, C..1.  McHugh.'  In  its  original  form  (Figure  1),  it  was  clearly  the 
(mechanically)  sijiiplest  lieat  engine  ever  invented.  When  heat  is  applied 
to  the  boiler,  the  water  column  in  the  tube  oscillates  backward  and  forward 
and  the  toy  boat  is  driven  forward  through  the  water. 

McHugh's  discovery  was  apparently  serendipitous,  but  in  a later  patent,^  he 
seemed  to  have  reached  a qualitative  understanding  of  the  phenomenon.  When 
the  water  column  liits  the  hot  boiler,  steam  is  made,  and  the  increase  in 
pressure  pushes  the  water  c.)lumn  down  the  tube  towards  the  exit.  The  pres- 
sure drops  as  the  steam  condenses  on  the  cool  wall  of  the  tube,  and  the  outv\'ard 
momentum  of  the  water  column  gives  time  for  the  steam  pressure  to  fall  well 
below  ambient  before  the  direction  of  motion  is  reversed.  So  the  water  column 
is  sucked  back  into  the  boiler,  and  the  cycle  repeats  itself. 

A net  propulsive  force  is  produced  because  the  flow  into  the  pipe  is  a "sink 
flow"  from  all  directions  in  the  2v  sterad  about  the  pipe  exit,  while  the 
flow  out  is  a unidirectional  propulsive  jet. 

A quantified  understanding  of  the  McHugh  phenomenon  is  desirable,  partly  be- 
cause it  enables  the  princijilc  to  be  applied  to  larger  scale  applications, 
such  as  simple  pumps,  boat  propulsors,  fluid  agitators,  and  oscillatory 
heat  pipes  (essentially  all  the  heat  supplied  to  the  boiler  is  discharged  in 
the  jet),  but  also  (for  example)  because  analogous  phenomena  ("chugging")  can 
cause  high  stress  levels  when  emergency  cooling  water  is  forced  into  an  atomic 
reactor.  Finally,  it  is  peripherally  related  to  the  dynamic  mixing  of  water 
and  hot  lava  or  metal  which  is  thought  to  be  the  cause  of  a number  of  cata- 
stropliic  explosions.^ 

Apart  from  McHugh's  second  patent,^  the  first  published  attempt  to  explain 
the  cycle  seems  to  have  been  made  by  Baker''  in  19.T2.  Baker  was  wrong  in  his 
explanation  of  the  cycle  (he  assumed  that  a bimetal  diaphragm  was  needed  in 
the  boiler),  but  correct  in  explaining  the  thrust  mechanism. 

Other  analyses  known  to  the  writer  are  by  Dickmann,^  Miller,^  and  MacKay.'' 

The  latter  shows  that  lioilcr  resiliency  is  not  essential  to  the  cycle.  MacKay 
says  in  part : 


"The  boiler  can  be  rigid  and  still  give  rise  to  these 
oscillations.  A nice  demonstration  can  be  had  with  a 
Pyrex  tube  having  a bulb  on  the  end  acting  as  the 
boiler.  Eitlicr  small-amplitude,  almost  sinusoidal 
oscillations  of  the  column  can  be  observed  or,  by 
changing  the  heating,  violent  cyclic  emptying  and 
filling  of  the  bulb.  In  this  latter  case  the  sudden 
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I roscncrativc  condensation  of  tlic  steam  is  observed 

most  clearly  since  there  is  no  flexible  top  to  act 
as  a cushion.  I'or  best  oper.ation  all  air  should 
first  be  driven  out  by  boiling.  Professor  I- rank 
Goyan,  who  joined  in  several  of  the  recent  experiments 
has  emphasized  the  necessity  for  liquid  water  to  con- 
tact the  hot  region  to  provide  a supply  of  steam  for 
later  condensation.  This  factor  may  acc;ount  for  the 
somewhat  more  stable  operation  of  a two-tube  unit 
over  one  with  a single  tube." 

The  observation  about  removing  gases  from  the  water  by  boiling  is  particularly 
interesting  to  the  writer  in  the  light  of  our  more  recent  experimental  work, 
described  later  in  this  paper.  Gas  released  by  boiling  in  the  boiler  can 
become  trapped  there,  and  inhibit  the  engine's  operation.* 

MacKay'^  (who  was  apparently  unaware  of  Dickmann's^  work)  also  gave  a correct, 
if  qualitative  explanation  of  wlty  the  engine  developed  thrust.  But  Dickmann's 
was  much  more  detailed,  and  included  experimental  data  (Figure  2)  from  an 
apparatus  in  which  a water  column  was  oscillated  by  a crank-driven  piston. 

The  tlieoretical  line  in  Figure  2 is  for  sinusoidal  piston  motion  and  is  given 
by 

I Average  thrust  = T^^  = i p(wX)^A  (1) 

where 

u is  the  crank  speed  in  radians/scc 
X is  the  piston  amplitude  or  half  stroke 
A is  the  cross-sectional  area  of  the  duct 
P is  the  water  density. 

This  very  simple  result  was  derived  independently  in  References  8 and  9 using 
entirely  different  methods  of  analysis.  Siekmann^  integrated  the  total  momen- 
tum flux  in  the  vortex  trail  - an  involved  calculation  - while  the  writer® 
considered  only  tlie  (one-dimensional)  flow  within  the  duct.  The  latter  approach 
seems  preferable  in  that  it  allows  skin  friction  and  inflow  loss  effects  to  be 
included.  The  fact  that  the  experimental  thrust  is  slightly  lower  than  the 
theoretical  result  may  be  due  to  a number  of  causes: 

(a)  Tlic  reference  8 theory  assumes  "sink"  flow  in,  "jet" 
flow  out.  Dickmann® suggests  that,  right  at  the  start 
of  the  exhaust  stroke,  the  flow  is  "source- like"  in 
nature,  and  becomes  "jet-like"  as  the  vorticity  builds 
up  at  the  boundary.  (His  two-dimensional  flow  analog, 
model  of  the  downstream  flow  is  given  in  Figure  3, 
together  with  that  of  Sickinann's  later  analysis  in 
Figure  4.)  It's  possible  that  some  measurable  loss  of 
rearward  moment  occurs  as  a result  of  such  transient 
".source- 1 ike"  flow. 


♦However,  Finnic  and  Curl'*®  say  of  the  Figure  1 configured  engine,  "...  the 
performance  of  the  boat  is  generally  improved  by  the  presence  of  small  amounts 
of  air  in  the  chamber."  The  reason  for  this  disagreement  is  currently  quite 
unclear. 
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AVERAGE  THRUST  IN  GRAMS 


CRANK  SPEED  IN  REVOLUTIONS  PER  MINUTE 


Figure  2.  nquation  (1)  compared  with  the  experimental  measurements  of 
Trunz,  as  reported  by  Dickmann.^ 
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Figure  3.  Dickmann's^  two-dimensional  flow 
model  of  a pulsejet  wake. 

(a)  Induction;  (b)  Discharge. 


Figure  4.  Sickmann's'*  two-dimensional  flow 
model  of  a pulsejet  wake. 


(bl  [Equation  (1)  assuiiios  that  tticro  is  no  inlot  pressure 
loss.  It's  not  clear  (lo  the  writer  at  least)  what 
tlic  effect  of  a flow  contraction  prc:ssiiro  loss  would  1)C 
during  an  instrokc.  Tlie  avorajjc  tlirust  would  he  less, 
because  the  intake  thrust  wouUl  l)c  ncju'itive  instead  of 
zero.  But  the  wake  is  jjenorated  by  the  outstrokc,  so 
we  should  tlien  have  the  anomalous  result  of  the  thrust 
being  less  than  the  momentum  in  the  wake. 

(c)  Equation  (1)  assumes  that  the  D'Arcy  friction  factor  f 
is  constant,  so  that  skin  friction  terms  cancel  out. 

The  magnitude  of  the  error  which  this  gives  has  not 
yet  been  tested,  but  one  would  not  expect  it  to  bo 
large,  since  inward  and  outward  velocity-time  histories 
arc  identical. 


Cavitation  does  not  seem  likely  in  the  Dickmann  experiments, 
pressure  would  be  about  150  Ib/ft^  below  ambient. 


since  the  minimum 


Some  further  discussion  of  lUckmann's  work  was  given  by  Schuster,  et  al’*^  in 
1960  (translated  1963),  still  from  the  vantage  point  of  determining  thrust  by 
integrating  the  momentum  in  the  wake.  In  1963,  Finnic  and  Curl  (who  were 
aiijiarently  unaware  of  Dickmann' s work)  presented  an  apparently  rigorous  analysis 
to  show  that,  in  irrotational  flow,  the  average  thrust  during  the  instrokc 
must  be  zero.  For  the  outflow,  they  obtain  the  conventional  result 

_2 

Thrust  = pAu 


that  is,  the  product  of  the  instantaneous  mass  flow  rate  and  the  mass-averaged 
velocity. 

A review  of  previous  work  would  not  be  complete  without  mentioning  that  of 
Gongwer'°>'’  at  the  Aerojet-General  Coi-poration , even  though  he  was  not  using 
the  McHugh  cycle.*  A typical  Aerojet  pulsejet  is  illustrated  in  Figure  5. 

Fuel  is  burnt  in  the  combustion  chamber,  and  the  expanding  gas  forces  the 
water  out  of  the  tailpipe  to  form  a propulsive  jet.  As  the  pressure  drops, 
a fresh  charge  of  water  is  admitted  by  the  one-way  valves  and  the  cycle  is 
repeated. 

Efficiency  was  found  to  be  fairly  low,  presumably  because  it  is  what  used  to 
be  called  a "Type  I"  engine^ ^ without  precompression.  This  roughly  halves 
the  cycle  efficiency.'^  Valve  life  was  short  under  the  repeated  mechanical 
loading,  and  vibration  was  very  high.  Figure  6 sliows  a typical  pressure-time 
history,'^  and  explains  the  latter  problem. 

Principally  for  these  reasons,  development  of  the  Aerojet  pulsejets  was  even- 
tually abandoned.  But  the  work  did  demonstrate  the  feasibility  of  developing 


* Of  course,  some  steam  must  have  been  generated  at  the  gas/water  interface. 
Heat  lost  to  the  water  would  reduce  efficiency,  but  that  used  to  make  steam 
would  not  be  a total  loss. 


6 


8 


thrust  forces  as  as  1000  lb  at  00  knots  forward  speed,  and  that  gas 

could  impose  large  acceler.it  ions  on  a much  heavier  fluid*  without  being  con- 
figured as  a Hunijiiirey  gas  purnp.'’’ 

The  writer  became  inten-sted  in  the  McHugh  cycle  in  1971,  and  we  have  since 
constructed  and  tested  a number  of  engines.  Iv'e  ]irefer  the  generic  name  water 
pulsejcts  for  any  pulsed  proimlsor'®  and  steam  water  pulsejets  (SV.'I’J)  for 
those  having  the  McHugh  cycle. 

F'igure  1 shows  our  basic  engine  conceptually.  Rather  than  being  larger  than 
the  pipe,  the  "boiler"  is  cither  the  same  size  or  smaller.  (A  typical 
laboratory  engine  has  a pipe  diameter  of  one  inch.)  The  boiler  wall  is 
thick  so  that  it  can  capacitively  discharge  heat  to  the  water  during  the 
few  milliseconds  that  it  is  in  the  boiler  cavity.  It  is  heated  by  electric 
resistance  elements  for  experimental  convenience,  although  propane  gas  flame 
heated  boilers  have  also  been  employed  successfully. 

Figure  8(.'0  shows  the  stcam/watcr  interface  travelling  away  from  the  boiler 
in  a transparent  duct,  and  the  opposite  motion  ii.  Figure  8(b).  Note  that 
a "stripped  boundary  layer"  is  left  behind  by  the  outgoing  water,  and  picked 
up  on  its  return.  Although  not  very  important  hydrodynamically,  this  stripped 
water  is  thought  to  dominate  the  thermodynamics  for  much  of  tlie  cycle. 

Figure  9 shows  a selection  of  boilers  that  have  been  evaluated  experimentally, 
the  boilcr/pipe  joint  being  shown  in  Figure  10.  The  assembled  engine,  shown 
in  Figure  11(a),  has  been  sectioned  to  show  the  internal  details,  and  may  he 
compared  with  the  "multi-hole"  boiler  drawing  in  Figure  9. 

In  addition  to  the  basic  thrust  engine,  we  have  built,  during  the  last  six 
years,  a number  of  alternative  configurations.  The  bifurcated  tailpipe  of 
Figure  11(b)  is  intended  to  permit  ram  recovery  when  tlic  engine  is  operating 
at  higli  forward  speed  in  the  water.  In  the  present  program,  ram  iiressures 
up  to  374  Ib/ft^  were  imposed  on  the  inlet  side  (equivalent  to  19.4  ft/sec 
with  lOOv  ram  recovery)  without  any  significant  change  in  engine  oiicration. 

A suction  pump,  shown  in  Figure  11(c),  is  an  obvious  extension  of  such  a 
bifurcated  thrust  engine.  To  date,  such  pumps  have  only  been  operated  in  a 
suction  mode'*^  but  since  a thrust  engine  works  against  high  back  pressure, 
when  a nozzle  is  fitted  (about  8000  Ib/ft^  for  a nozzle  area  ratio  of  three, 
corresponding  to  a lift  of  130  ft),  there  seems  to  he  no  fundamcnt.nl  reason 
why  they  should  not  pump  against  a head. 

I.uberoff**^  has  pointed  out  that  the  basic  (nonhi  furcated)  pulscjct  is  an  effi- 
cient agitator  which  may  have  applications  in  the  chemical  industry,  parti- 


* It's  sometimes  thought'^  that  Taylor  instahi lity'^* precludes  accelerating 
a fluid  with  a gas,  but  Taylor^^  and  l.ewis*^  did  not  say  this.  The  insta- 

bility results  in  fluid  "tentacles"  reaching  out  into  the  gas,  but  as  Lewis 
observes  "in  spite  of  these  very  large  surface  disturbances,  the  main  body 
of  liquid  ...  is  accelerated  as  though  they  did  not  exist. 


’bOILER" 


WATER  LEFT  BEHIND  BY 
THE  INTERFACE  ^ 

(a)  Discharge  Stroke.  (For  inviscicl  flov;,  and  no  nozzle,  the  equation 


of  motion  is  p -p 


UP  BY  THE  RETURNING  INTERFACE 

(b)  Induction  Stroke.  (For  inviscid  flow,  the  equation  of  motion  is 


Figure  7.  The  simplest  exhaust  and  intake  flows,  (p  = mass  density  of  the 
fluid.) 
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Fijiurc  8(;0.  On  the  outstrokc,  the  steam/wntcr  interface  is 
moving  from  left  to  right,  leaving  a "stripped 
boundary  layer"  behind.  On  this  engine,  the 
internal  diameter  of  the  le.xan  pipe  if  5/4  in. 


Figure  8(b).  Moving  right  to  left  on  the  instrokc,  the  interface 
is  picking  up  the  stripped  boundary  layer  water  as 
it  goes.  The  latter  has  Iiad  time  to  more  or  less 
fall  to  the  bottom  of  tlie  tube. 
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Frames  1-4 


Frames  5-8 


Figure  8(c).  In  this  composite  from  a high  speed  movie  (roughly  1000  frames 
per  second),  we  have  printed  every  fortieth  frame.  In  the  first  one,  the 
steam/water  interface  is  just  about  to  emerge  into  the  transparent  section. 

In  the  second  frame,  the  lower  half  of  the  tube  to  the  left  of  the  interface 
appears  half  full  of  water,  but  this  is  an  illusion;  the  "stripped  boundary 
layer"  is  merely  coating  the  inner  wall,  and  the  inner  surface  is  sufficiently 
disturbed  for  it  to  appear  opaque.  In  the  third,  fourth  and  fifth  frames, 
we  see  the  boundary  layer  water  fall  to  and  puddle  in  the  bottom  of  the  tube, 
while  the  interface  motion  is  arrested.  In  frames  six  to  eight,  the  interface 
moves  back  to  the  boiler,  picking  up  the  puddled  boundary  layer  water  as  it 
goes.  In  the  movie,  the  vorticity  generated  by  this  is  clearly  shown  by  the 
motion  of  the  gas  bubbles  behind  the  interface. 


STRAIGHT  ANNULAR  BOILER  MULTI-HOLE  BOILER  OPEN  HOLE  BOILER 

0.1  INCH  ANNULUS 


ELECTRIC  HEATER  CAVITIES 
(SAME  ON  ALL  BOILERS) 


COPPER  CENTER -BODY 


Figure  9(a).  Boiler  configurations  tested. 
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Figure  9(b).  Boiler  configurations  tested. 
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Figure  11(c).  A pulsejet  punp.  As  for  the  Figure  11(b)  configuration,  one-way  valves  are  fitted 
to  the  foot  and  discharge  ends.  Suction  lifts  of  up  to  11-fcct  have  been  achieved 
without  significant  diminution  of  the  flow  rate. 


cularly  where  corrosive  fluids  mitigate  against  the  use  of  moving  parts,  seals, 
etc . 

In  what  follows,  wo  shall  confine  ourselves  to  the  simple  thrust  engine,  al- 
though the  tlieory  will  generally  be  api>licable  to  a bifurcated  one  as  well, 
so  long  as  ttie  legs  are  of  the  same  length. 
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i:;vi.\r:axs  or  motion  for  an  iwiscio  i-i.tiiD  column 


The  Rasic  Fqiu^tior^ 

I'or  siniplicity,  v.'c  first  consider  the  inviscid  flow  case.  Then  if  p is  the 
static  pressure  at  the  piston  (or  steam  water  interface)  and  p^,  is  tlic  pres- 
sure at  the  noi'.zle,  Newton's  law  gives 

A(p.^  ■ ^ = mx  + mx  (2; 

where  A is  the  cross  sectional  area  of  the  duct. 

Now  the  fluid  mass  is  m = pAx 


m = pA.x 


for  the  instroke 


For  the  outstrokc,  m = 0 because  momentum  is  conserved  in  the  jet.  So  we 
have,  for  the  inviscid  flow  ecpiations  of  motion 


Instroke 


Pn  - P 


= XX  + (x)^ 


Outstrokc 


Px-  - P 


During  the  outstrokc,  p^^  = p^,  the  ambient  value.  The  same  is  true  during  the 
instrokc,  but  the  reasoning  is  more  complex  and  lacks  rigor.  From  Bernoulli 


1 34) 

P»  - P aT 


where  is  the  velocity  potential.  Now,  as  will  be  shown  below,  equation  (3) 
gives  zero  average  force  on  the  instroke  for  = p^,  and  this  is  clearly  the 
correct  result,  since  all  fluid  is  stationary  at  the  start  and  the  finish  of 
a single  inward  stroke.  So  from  equation  (5) 

at  2 ^ 


(dx/dt)  dt 


i/" 


which  is  a valid  solution. 
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Ihc  proof  that  tho  avorap.o  instrokf  thrust  is  zero  is  simple,  since  from  (3) 


‘.t  ,..T 

1 N n \ 

Ap 


:/■ 

IN 


(p  - p)  cit  = 


/ XX  dt  + / (x 

IN 


. ? 

x)  dt 


Hvaluatint;  the  first  integral  by  parts  gives 

■At  ,.At 


At  ...T 
IN  ay 

Ap 


XX 


f .2  r.2 

- / (x)^  dt  + / (X)^ 


dt 


since  X = 
tlie  same 


0 at  the  beginning  and  end  of  the  induction  stroke, 
analysis  for  the  oiitstroke  (ocpiation  4)  gives 


t, 


t T 
P av 


I- 


pA  / (x)"  dt 
At 


Not  ice  that 


which  simply  says  that,  for  constant  duct  area  A,  the  force  impulse  (tpT  ) on 
the  end  of  the  boiler  is  equal  to  the  time  integral  of  the  momentum  disefiarged 
in  the  jet.  This  is  the  result  obtained  by  I’innie  and  Curl.''® 

(In  the  above  equations.  At  is  the  total  instroke  time,  and  t.,  the  total  cycle 
(in-  plus  outstrokc)  time.] 

The  same  type  of  analysis  can  be  used  to  show  that  the  work  done  is 


4 


dt 


where  the  plus  sign  refers  to  the  instroke  and  the  negative  sign  to  the  out- 
stroke,  for  which  x < 0. 

When  the  nozzle  area  (Aj,j)  is  not  the  same  as  A,  for  the  oulstroke,  we  h!»ve, 
from  Ilernoulli* 

Pn  = I"™ 

where  . 


*In  this  case,  the  3(t)/at  term  is  the  inertial  component  which  is  already  inclu- 
ded in  the  equation  of  motion. 
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(7) 


So  t '10  i lake  the  f'ina>  form 

!\„  - I' 

- - XX  + a(x)2  |>  (say) 


wiicro 

u = 1 for  tlio  insnaike* 

“ -r,  for  the  oiitstroko 

If  P is  constant,  we  can  solve  this  in  the  convent i onn I phase  plane  by  writing 
u = X,  X = 

dx 


I* 


du 

xii  + au 
dx 


Then 


■ / 


2 

^2 


dx 

X 


"1  ''1 
which,  after  some  reduction,  becomes 


'^2  7i  - (Xj/X2)^‘*]  + (Xj/X2)^‘*  (a  0) 


= 2P  log  (X2/Xj)  + Uj  (a  r.  0) 


(8) 


(9) 


A second  integration  is  required  to  find  x = f(t)  and  this  is  analytically 
])ossiblc  only  for  certain  values  of  (a);  namely 


o - ,1111 

® ^ ’ 2 ’ 3 ’ ^ ’ 5 ’ ‘ ' ’ 


These  solutions  can  be  obtained  by  the  transformation 
z*"  = [(P/a)  - uj  ](xj/x2)^“ 

whore  m = 1 or  2. 

A different  method  of  solving  equation  (7)  involves  the  ti'ans format  ion 

n 

X - z 

* The  same  result  as  Finnic  and  Curl***’  but  obtained  in  a very  different  way. 
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so  th;it 


X = n 


n-1 


, ,,  n-2,-.2  . n-1. 

n(n-I)z  (z)  + nz 


Makiiiji  tlicsc  .suhst  i tut  ions  in  equation  (7) 

, 2(n-n,.,2  2n-l..  2 2(n-l),.,2 

n(n-l)z  (zj  + nz  z + an  z (z) 


= P 


Tile  (z)”  terms  will  cancel  out  if  n = l/(l+a),  giving 
1 fl-a)/(l+a) 


Ua  ^ 


z = P 


(10) 


Writing  q = z,  ~ ^ *1 


f"'  r' 

J q dq  = (1+a)  / P dz  \ 

\ ( 


1,2  2 , 
2 ^^2  - '’l  ^ 


(11) 


The  inverse  transformation  is  given  by 

1/n  a+1 

z = X = X 


d z , , ^ a 


X 1/ (n- 1)_  . a 

--  z - (a+llx  X 

n 


,2  2a  , . 2 2a 

(x^)  x^  - (Xj)  Xj  = 2 


2 

/r 


dx 


‘1 


(12) 


If  P is  not  constant,  it  often  varies  as  a function  of  x;  as  in  isothermal 
expansion,  for  example.  So  if 

2a- 1 


dx 


is  intcgratibl e , we  have  a solution  for  (x)"  in  terms  of  x.  This  is  often 
sufficient  for  performance  calculations,  where  is  zero  at  the  end  of  the  oiit- 

X, 


St  roke , i . e. 


, ,2  2a  • 

dx  + (..j)  Xj  =0 


(13) 
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1 .)!•  t !k' 


S' ire  case,  tliis  leatls  to 


r , :a  2a  , ..  ^2  2a 


^2  . > 1 / '^a 

. [ 1 - afx^  )“/l*l  " 


(14) 


This  c'C|iiitii)n  is  tiie  l)asis  of  the  sinijilo  piilscjot  cycle  theory  presented  later 
in  the  paper. 


The  "Hr ink  ini;  Str.aw"  Pruhlern 

faMiiliav  a])plication  of  tliesc  equations  is  the  notion  of  the  liquid  in  a 
di'inlvinj;  straw.  If  the  tipper  end  is  held  closed  wliilc  it  is  thrust  vertically 
into  ;i  liquid,  tite  fluid  inside  the  straw  will  be  depre.s.scd,  as  shown  in 
fiyure  12.  Ileiiiovin!'.  tlie  obstruction  at  the  upitc-r  end  will  .allow  tlic  fluid 
to  rise,  and  in  this  case,  the  equation  of  motion  is 

XX  + = g(h  - X)  (If,) 

Substituting  P’  = g(h-x)  in  equation  (13)  and  noting  that  a = 1 for  an  "instroke" 
and  Xj  - 0,  wo  obtain  the  following  relationship  for 

-X 

2g  / (hx-x^)  dx  = 2g  f I (x,^  - J (^m  • ° 

Writing  x = ~ h/Xj 

x^  + (1  - I h)x  + (1  - I h)  = 0 (16) 

This  result  is  plotted  in  Figure  13.  IVlicn  tlic  top  of  the  straw  is  released, 
the  fluid  in  the  straw  will  rise  above  the  free  surface  to  the  height  given 
in  Figure  13,  before  falling  hack.  (Actually,  somewhat  less,  because  of  vis- 
cosity, of  course.) 

The  downward  motion  will  he  governed  l)y 

XX  = g(h  - x)  (17) 

for  whicli  the  solution  is 

h log  X + (1  - x)  = 0 (18) 
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0 2.0  4.0  6.0  8.0  10.0 

TUBE  IMMERSION  PARAMETER  h/x, 

FiKiire  1 .V  The  overshoot  of  .-.i.  initially  depressed  column  of  inviscid 
fluid  in  :i  drinkiiiR  straw,  from  ecpiation  (l(')- 


As  h'.'fori’  li 


initial  x 


siinimuiii  X 
initial  X 

The  solution  to  cc{iiation  (18)  is  jilottoj  in  Fij^iiro  14.  As  an  example,  suppose 
a straw  were  iiniicrsed  two  inches,  with  the  inviscid  li(|uLd  inside  at  a height 
of  two  inches  above  the  free  .surface  (h  = O-.h).  Then  relative  to  the  free 
surface,  the  successive  maxima  and  minima  would  be 

+ 2.0  inches  (the  start) 


+ 0.8.T  " 

- 0.04  " 

+ 0..A1  " 

This  pre.scnts,  thercT'ore,  the  rather  rare  spectacle  of  damped  motion  of  an 
inviscid  fluid. 


Sinusoidal  Motion  Forced  by  a Piston 

Equation  (7)  is  readily  solved  if  the  motion  is  prescribed.  Let 


X - Xq  + X sin  0 


6 = «t 


so  that 


Also,  let 


X = loX  cos  0 


X = -(jj  X sin  9 


Y - XyX 


Then,  for  the  inst  rokc 


P - l’„ 

'P  ” T 2' 
IN  pm  X 


2 2 

Y sin  0 + sin  0 - cos  0 


Tl  3 

For  the  out  stroke  ( -^  < 0 < —it  ) 


2 2 

Y sin  0 + sin  0 + 0 
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rij;iir',s  1 fi  ana  lb  prescjiit  ('quations  (1'))  and  (20)  as  a fiiiu-tic'n  of  phase 
0,  and  as  eo:'..a  -itional  pressiti'o  volinnc  (P-V)  diagiam.  The  average  pressure 
on  the  piston  is 


instroke 


0 


(21 


out  stroke 


{(1  . o 


( 2 a 


To  euinjvate  the  thrust,  we  compute  the  impulse  (I)  of  the  discharge.  By  con- 
tinuity, the  jet  velocity  is 


p(A^/A^) 


Sti/2 

f : 

I cos 

kn 


cos  0 do 


uX 


ii 

1 


Since  the 

period  tp 

= 2ti/u 

, the  average 

thrust  corresponding  to  this 

i mini  1 sc 

is 

T = 

1 . 

J 

1 pA^(<..X)^ 

4 A“-  - 

c = 

_1n__ 

= J (A/A^) 

= {vrrr2c 

(23 

poj^X^A 

Comparing 

(22)  and 

(23) , we  see  that  the 

average  tlirust  force  on  t !ie 

no 7.3  1c  i 

AC 

,\02 

1 , 
'4 ' 

v^l  * 2c  - (J 

» O] 

(24 

A('|,  0 for  ' 0,  and  is  otherwise  negative,  as  would  be 

.\0Z  e.xpectcd. 


PISTON  PRESSURE  PARAMETER  Ap//>0) 


PRESSURE  COEFFICIENT 


_ VOLUME 

VOLUME  COEFFICIENT  Cw  = 5 

ttR  2x 


Figure  16.  liiiuMtions  (19)  and  (20)  proscnted  as  a "pressiirc-voli'me" 

diagram.  The  area  between  the  curves  is  proportional  to  the 
work  required  to  drive  the  piston  through  one  cycle.  (c=0) 
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riu'  work  iwwji' i ’Wwl  to  move  t lio  piston  is 


(:ompririn}>  (26)  ;iritl  (27),  wc  see  that  tlic  work  done  hy  the  piston  to  expel  the 
fharge  appears  as  kinetic  energy  in  the  jet,  thus  validating  our  etpiation  of 
motion  foi-  the  discharge.  Ttie  work  done  by  the  piston  to  fill  the  duct  does 
not  appear  as  energy  in  the  jet. 

Ihe  result'^  can  readily  lie  extended  to  the  case  where  the  duct  is  moving  through 
a fluid,  the  end  heing  hifurctited  to  recover  ram  )>rossure,  or  to  pump  from  one 
pressure  to  another.  It's  found^'^  th.it  the  ideal  propulsive  efficiency  is 
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K)icrt' 


: (1^0 


2'. 


1 

"> 


1.  ^ 11  /ioX 

o 


u fret'  strernn  velocitv 

o 


f28) 


I'or  r,  - 0,  Ko  }',ct  tlic  low  v'aluc  of  - 0.267  at  u = 0.4.S5.  I:lTicieiicy  in- 

crcasc.s  with  c to  a limit  of  0.7266  as  C -*  1'h  i s strange  result  c.aused 

some  pujirlcmcnt  until  it  was  rc.ilizcil  that,  for  part  of  the  sinusoidal  strol.e, 
the  jet  velocit>'  is  less  than  that  of  the  fveestroani.  Increasing  the  nozzle 
cnnt  ract  i on  ratio  (increasing,  r)  reduces  tlie  fraction  of  the  stroke  when  this 
is  so,  and  hence  improves  efficiency. 


'Ihis  led  to  an  analysis  of  the  effect  of  jet  velocity  nonuniformity  on  propul- 
sive efficiency,  reported  in  Reference  ,8  from  whicli  figure  17  is  reproduced. 

It  was  found  that  a uniform  jet  velocity  results  in  the  same  propulsive  effi- 
ciency as  a I'roudc  actuatcir  disc.  The  more  the  velocity-time  liistory  departs 
from  this,  the  lower  tlu'  j'ropulsive  efficiency. 


The  I mjm  Ij;  i.v<L J’"  lit? 

(•igurc  IS  shows  the  ty[)e  of  steam  pressurc-t  ime  history  measured  with  a steam 
water  iMilsejot;  and  of  course,  the  tlirust-tiiiie  history  looks  very  similar. 

A simple  idealization  of  this  is  to  assume  that  the  water  column  motion  is 
impulsively  reversed  at  the  boiler,  and  that  at  all  other  times,  the  pressure 
is  constant  at  some  value  p^.  .So,  apart  from  the  impulse,  the  motion  of  the 
water  is  described  by  equation  (7). 


Let 


Then  , 


a 

from 


bo  the  boiler  ordinate 

the  point  at  which  the  interface  is  farthest  from  the  boiler 
the  water  velocity  jirior  to  boiler  impact 
the  water  velocity  immediately  after  boiler  impact 
equation  (9),  for  the  foke  (a  = 1,  iij  = 0) 


Till  - (.„o„)h 

[>’  = (Po,  - P^A'I 


(29) 


32 


STEAM 

PRESSURE 

Ps 


AMBIENT 

PRESSURE 

Pcc. 


VACUUM 


Fipiiro  18.  Idealized  and  observed  steam  pressure  in  a pulsejet  (not  to  scale) 
as  a function  of  time. 
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i'C>r  the'  (Mit  St  rc'ko  fn 


x.j),  Vc  havo  from  fvtuation  (11) 


X 


X 


B 


. ? - I / '^r 

n ^ ^(x.^)Vi>i 


{r.  i 0)  (30] 

(■•.  - 0]  (31) 


'CIr'  Cl' ! i\-sj)oiul  i nr,  Imilor  exit  velocities  are 


-yii'A)i(X|,/x„)-''  - 1 1 


/2i'  lot  (x„/ri 


(c  if  0) 
(r.  = 0) 


(32) 

(33) 


So,  knowiuj;  and  Xj^,  ciiuations  (20),  (32)  and  (33)  tell  us  i;hnt  the  iminilsivc 
velocitN'  change  (x  - x,^)  at  llu'  boiler  must  be  for  a given  value  of  P.  The 
inijHilse  assoc  i atod  wi  th‘ th  i s is 


mAu  = 

pAxB(ip  - 

f)AX|ji'T^ 

■ 'V'li’'  * 

/n/c)[(xg/.xj‘-''’  - 11  j 

j 

(C!^0) 

(31) 

pAXg>^ 

y/l  log  (x^/.x^)  j 

(r.=0) 

(33) 

During  the  rest  of  the  cycle  (in  this  model),  tlic  pressure  is  below  ambient, 
so  that  the  total  iminilsc  per  cycle  will  be  less  than  this  value,  be  now 
determine  the  cycle  time.  For  the  ijis^rol^',  we  have  from  equation  (9)  (u^  = o) 


i = vA’ll  - (x„/x)^l 


)" 


(30) 
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.ir-  t!i;'  O'lt V.m’.o  . I'l'oin  (.iiiiat  ion  (12)  (;i  - - , x,  - x , x,  =-  x„) 

— I 4I  1 n 


(xlVx'"  - 


- 7 r 2 r 

-(P/0|x  - x,j 


r . 7 2r. 

■ 1-0:1  L'i|u;irio?i  loJ'  - (x,  )“  = ‘ * 


{.  ,-,2  -2r. 
,,  (xl  X 


-2r,  -2r. 

X - X 

in 


'/  (!'/.' 


(!Vc)[(x/x^)-'  -1] 


r . - 

J /(x/x  - 1 


- A’/t.  I dr 


At  y'P/r. 


Analytical  Holutions  arc  po;;siblc  for  A = 1/-1,  1/2  and  1,  and  arc  given  in 
Talilc  1 and  I'ij’.urc  10.  N'uncrical  integration  for  other  values  of  C arc  given 
in  Table  2.  The  Form  of  Tnhlc  1 is 


, , A' 

T,  ~ ^OUt  X., 


X /•' 

-J"-  / 

Hi 


/^x/.xj^^  - 1 


x.,/x 
B m 


Tor  the  case  of  c,  = 0,  we  have,  from  equation  (9)  (Uj  = , Xj  = x^^) 


-/2V  /og 


log  (x/Xj,)  + 


But,  from  equation  (S’^) 


(’‘a^  _ , ’'b 

2P  ■ X 

m 


X - - 


21’  /log  x/x 


-C 

I /leg 


= /2P  At 
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Tahlo  1.  Analytical  Solutions  to  liquation  fSS)  for  1^ 


1 o;; 


, 2"  , 

'B'  -'Ki*  ■ * X 

ni 


] 


. 2 -“T 

B m 


1 

r. 


2/2  -11-^  (x„/xj^^^ 


B m 
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T,iM. 


. 1 
.? 
.3 

.5 

.e 

.7 

.p 

,o 


.599  3 5f 
.79P1?7 
.9)6^179 
.9P513 
1 .(?07/J 
.97  6P7P 
.91-!i535 
.81391^ 
.61P513 


r> 

<• 

? 

S' 

? 

p 

p 

p 

p 

p 


4 

4 

4 

4 

4 

4 

4 

4 

4 


. 1 
.P 
.3 
.4 
.5 
.6 
.7 
.P 
.9 


.4p?6f 

.646513 

.7949  18 

.R916P3 

.938577 

.928437 

.884559 

.79RP75 

.6  07  {'58 


.1 

.P 

.3 

.4 

.5 

.6 

.7 

.8 

.9 


.31046 

.5294PP 

.69 1539 

.807398 

.874301 

.8828 88 

.855368 

.782464 

.601635 


. 1 
.2 
.3 
.4 
.5 
.6 
.7 
.8 
.9 


.241565 

.442502 

.606345 

.733215 

.815084 

.839507 

.827051 

.767102 

.596247 


5 

5 

5 

5 

5 

5 

5 

5 

5 


. 1 
.2 
.3 
.4 
.5 
.6 
.7 
.8 
.9 


. 198329 
.378795 
.537397 
.668892 
.761152 
.798647 
.799685 
.752007 
.590894 
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aM.' 

1 aoi!  t ! !iu'a.!  I 

) - 

x 

1, 

R ll 

c 

6 

. 1 

.169880 

6 

.331765 

6 

.3 

.481074 

6 

.A 

.613681 

6 

.5 

.7  1 P47  1 

6 

.6 

.76P352 

6 

.7 

.77  333 

6 

.8 

.737 196 

6 

.9 

.58558 

7 

. 1 

. 1 5PP8  1 

7 

.? 

.P964  19 

7 

.3 

.437354 

7 

.566534 

7 

.5 

.6688 16 

7 

.6 

.7P4648 

7 

.7 

.748829 

7 

.8 

.7PP684 

7 

.9 

.58P3P4 

8 

. 1 

.13559? 

8 

.P 

.P69P37 

8 

.3 

.4P1 175 

« 

.4 

.526313 

8 

.5 

.629831 

8 

.8 

.6915G4 

8 

.7 

.723815 

8 

• 8 

.708486 

8 

.9 

.57  50  69 

9 

. 1 

.124501 

9 

.S 

.247816 

9 

.3 

.37153? 

9 

.4 

.491933 

9 

.5 

.59  5003 

9 

• 6 

.660837 

9 

.7 

.700699 

9 

.8 

.694613 

9 

.9 

.56987  5 

1 fi 

. 1 

. 1 1 5699 

18 

.P 

.230533 

18 

.3 

.346947 

18 

.4 

.4624  24 

18 

.5 

.564156 

18 

• 6 

.632532 

IB 

.7 

.678683 

IP 

.8 

.68 1074 

IP 

.9 

.564725 
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Aj;ain,  thir.  must  he  intcgraf  cd  miinuiical  ly,  and  this  !ias  been  done  in  '/'able  3, 
in  the  format 


I 

o 


'■''="-.>'''011  r 


d(x/x  ) 

_ fi. 

/loK  (x/x^l 
in 


l-'roin  eciuatior.s  (36),  (3SJ  and  (40i  , the  total  period  is  yiven  by 


(40) 


2^.  t 


/T- 


(c  i 0) 


2 . 1 111 


m-  'b'  W2  r,“ 


(C  = 0) 


(41) 


F.xccpt  for  tlic  case  of  r;  = 0,  we  cannot  compute  the  thru'-t  from  internal  pres- 
sures, because  the  force  reacted  hy  the  nozzle  is  (at  jircscnt)  unknown.  It's 
therefore  necessary  to  compute  the  momentum  in  the  jet.  Repeating  the  rea- 
soning employed  to  derive  et/uations  (23)  and  (35),  the  jet  impulse  is 


I.  = P(A^/A^,)  J (i)^  dt  = - P(A^/Aj^)  J i d.x 


= -pf-  X /pTc 
" n> 


= - pAx  /2P 

m 


r.  ^ /(x/x 

D m 

r\ 

I flOO  (x/x^ 

’'x„/ X 

B m 


(6  ^ 0) 


(f.  = 0) 


(42) 


For  the  case  6 = 0,  integration  by  parts  leads  to 


I . 
J 


pA.Xr  •'log  (x/x^)  - i I J 


(C  = 0)  (43) 


where  f is  defined  by  equation  (40)  and  has  already  been  tabulated.  Average 
thrust  is  then  given  by 


* Care  has  to  be  exercised  when  integrating  by  parts,  because  of  the  indeter- 
minancy  in  signs  between  l/x  and  log  x.  While  there  are  more  elegant  ways 
of  selecting  the  correct  sign  rather  than  +lil^) , a quick  check  on  a 

modern  electronic  calculator  is  perfectly  rigorous. 
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;iv  1. 


/_  (x_„/x,.)  I 
>'2  m B o 


lor  this  s;imo  case  of  f - ()^  wo  can  compute  the  thrust  from  the  inipulse  (equation 
as)  anti  the  internal  pressure,  i.e. 

T.. „ = r-  - pAP 


’f"l7V 

^ /2  ^\/S^  'o 


v.’hich  reduces  to  ctiuation  (It) 


For  the  case  c,  0,  wc  cannot  reduce  equation  (42)  to  the  form  of  I in  (38), 
so  we  have  to  solve  the  integral  nuinerical ly.  This  has  been  done  in  Figure 
20,  in  the  format 


x„/x 

B III 


T = pAx;^l>  = '^20/r: 

m 


x/x-J-'-l  d(x/xj 


so  that  the  average  thrust  coefficient  is  given  by 


C ^ 

(Xg/xJ[/i  - (x^/Xg)2  . i^J 


(Note  that  two  analytical  results  for  l.p  are 
1 t _ 4 ,3/2 


for  C = y 


3 


forr  ^ 1 I = f |(Xp/xJ  /(x,/x„)^  - 1 - log  / . 

( m ) 

One  could  now  proceed  to  compute  the  pressure  forces  on  the  closed  end  of  the 
engine,  as  was  done  before  ecpiation  (45),  and  by  subtracting  equation  (47) 
from  this,  determine  the  average  internal  force  on  the  nozzle. 
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I ft' Li' i on.,  y in  t tu-  I injui  1 s i vc  ('yo  k- 


.■\t  the  inst.it’.t  of  l)(>i  Ion  inip.u  t , the  innss  of  Kritoi’  in  the  thiel  is 
ener;;>-  .uldeJ  to  t li  i s mn-.s  of  when  it  k-rivrs  the  itoilei'  (with 

veloeity  .x.^l  is  therefore 


IV 


\ c.AX|,| 

2 ■ " ' <W"i} 

(r.o) 

;nA.y.]2  1,,,  (.yx^y  - |i  - (x„/.x„ih  1 (oo, 


I ro;;i  e([u;itions  (2hl,  (o-’i  ami  (e.A).  It  can  ho  showii^^  usinjt  these  cq'iat  ittits 
that  IV  is  oqii.'il  to  t Ix'  kinetic  cnerf,y  d i sch.'i  rj;eil  in  tlie  jet.  The  r;ost  ine.aniny, 
fnl  rie.asiue  of  "boiler  efficiency"  in  this  model  is  therefore 


"L 

"ll  * A(j.I 

'vhere  .'(j  i s the  heat  released  per  pulse 


(r,n 


I is  the  mech.'in  i c;i  1 equivalent  of  heat. 


So.  kno’.sinj',  X|||  .and  I’,  and  the  hoiit  input  to  the  boiler,  the  efficiency  c;in  be 

determined  tit  once  from  equations  (48).  It  will  be  an  inuler  est imate,  because 

to  .achieve  .a  y i von  stroke  (.x,,-x  ) with  a real  viscous  fluid,  the  boiler  exit 
, ^ B m 

velocity  x,  must  be  eroiitcr  than  needed  by  an  inviscid  fluid. 

u 


The  most  convenient  ptirameters  to  measure  arc  x and  the  period  tj,.  So  using 
equ.'itions  (41j  to  elimintiTc  I’  in  (SO),  wc  get 


2 - 11  - . qj"  (On)'j 


f'Ax, 


2 


!•-  ’"i’  ^ /a  .xf:  ’o( 


f 


B "I  > 

(C-O) 


In  the  ;;,anie  w.ay,  one  c.in  employ  equations  (41)  to  eliminate  P from  the  thrust 
equations  (44)  .and  (47),  giving 


ri  , ^ ^ 

" 1.,  1 *^i  ' (x  /x„r  ^ I 

^ ^ X|^  1 ' m It 


|v''l();.;  {y-J^  ) 
B i:i 


^ lit  t It  t ‘ I 

, ll>l  - (X  /xj  + u - 1 

y Xj^  o III  1!  /.?  X|^  o 


li^nrt'  21  shows  that  this  coirolnlos  cx|H'r  imoiit  a 1 data  c|iiitij  woll.  rho  l('ii<;ci’ 
Iho  stroke,  t ho  further  above  the  tlieory  I ini-  ai'c  the  expor  i nicnt  a 1 points,  .as 
would  bo  oxpooted  bcoaoso  of  tlio  j'.i'catcr  skin  friction  in  tiu'  loiij'.or,  fa.-'ter 
St  roko . 
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4.1 


CONSTANT  PRESSURE  INVISCID 
FLOW  THEORY  (EO.  51) 


Sdiiic  c I :i  I ilil.i  v.nni|M  fi'il  uilh  t i i'll  (^'). 

\s  ■ ■ \ii  1 ,1 1 II' 'll  ill  till'  tf.vl,  1 III-  il  i ;;i.l'i'|';iiK'v  liv-iui-i'ii  tlu'i'i'.'' 

• III'!  r \ |U' r i rii 'll  I I ■■  III' i IK' ! |i;i  I I iliu'  to  v i 1 1 v . (IIk'i>i'>  ,iiuI 
t l•■■l  il.it  :i  for  i'  0)  . 


I i 


Ri.Ai.  i i.iiiii  i-.nrcrs 


T!k'  "St  !■  i i'p;- J"  no!ii'<.I.'t7'>’ 

As  shtn.n  C'arii(.'i'  in  Fi};ur<.-  H,  tlic  slower  laoviiiK  tioiiiuln  ry  Inyi-r  water  is  "left 
hehiiu!"  on  the  ontstroke,  and  on  a long  ent'.ine,  has  time  to  collect  in  the 
hottoi’i  of  tliL'  duct  before  bein';  (licked  tip  by  the  returning  water  coliii:in.  A.s 
would  he  expectetl,  cons  i tleralil  e larp.e  scale  v'orticity  is  ileveloped  in  the  liead 
oi'  tlie  retiirnin;;  w.iter  column,  as  it  "rolls  over"  the  puddleil  boundary  layer 
wat  e r. 


l\'e  can  compute  the  (piantity  of  water  left  behind  liy  assuninp,  that  tiie  local 
veU'cit)'  in  tlie  (lipe  is  given  by 


(1  - r/lt) 


1/n 

•v  ‘ 


(SI) 


wbt're  r is  the  locttl  rtidius 

K is  the  internal  rtidius  of  the  duct 
u is  the  local  velocity 
V is  the  centre  velocity 

n is  an  ex|>onent  determined  from  experiment. 
Sch  1 ielti  ing-^ ' presents  datti  wliicit  imply 

a 1 I, 

n ~ 2.1  lie 


(5S) 


2n^V 


where  the  Reynolds  number  F’e  = 2uR/v 

u is  the  average  velocity  in  the  pipe  = (T+atyTl^^n) 
V is  the  kinein.it  ic  viscosity 

Relative  to  tlie  central  velocity  V e<iuation  (S4)  implies  a volume  flow  defect 
R -1 


'I 


1/n, 


.?T(r(V-uj  dr  ---  2ir\iR^  j (I  ■ {1  ••  r/R)  ' ](r/R)  d(r/R) 
o o 


2 

II RV 


(.SnH) 
(n^l j(2n+iT 


(B6) 


^ripjied  water  volume/unit  length 
itR^ 
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1'  1^'  i ’’  cor'.p'it  ,\1  'I'Mi;;  Iho  "lioilt-r  exit"  vi-locity  x , anj  assm^iing  a 20t)°l' 
K.itcr  t <- iiipo  ;■  1 1 .1  r-L'  iH-a  r tlic  i n i r‘rr,n.'>' , lluai  for  a typical  one-inch  piilsojct, 

V ;>.S  \ 10  ft  '/Nce,  he  ti  ■ 8.S4,  "SIVK"  “ O.lf., 

To  chock  this  result  o xjk  r i mont  a 1 ly  . a one-  inch  O.I).,  0.875  inch  1.1).,  pulse- 
jet  lois  built  uj)  froi:!  "1  I'.xpn''  tube,  iisinj’  the  tapered  circular  boiler  with 
cf  nter-heatcr  (Drc.  1 .1.5C:.>  I (.)  of  Piipire  '»(bj.  Total  duct  Icnp.th  wa.s  8.5. 2.5 
inches,  and  a sinijle  action  lecer  shut  off  valve  was  installed  .54.25  inches 
frc'ia  t hi'  boiler  attachment  flaupe  face. 

Uhen  the  euy.ine  was  ruuniiif'  stc.adilv,  with  all  temiK-vatures  stabilised,  niea- 
siirei'ients  of  thrust,  froipicsicy  aiul  boiK-r  tom[)eratiirc  were  made.  Tbrir.t  was 
Inca su |■■.'d  by  a force  cell  at  the  boiler,  and  also  deduced  from  exhaust  inomcntu:!, 
the  Vi-locity  boiiip.  mca-uired  with  a pitot  static  in  the  duct,  connected  to  a 
h i i'.h  Ireipjeniy  strain  yauy.ed  d i aph  rap.pi  d i f i'e.rent  i a 1 pressure  transducer, 
both  sipn.ils  were  strip  cliii't  rccordcil  and  siihseqiient  1 y integrated  by  h.inj. 

Ih,'  jet  i'clocity  result  was  ri-e.arded  as  tlie  most  accurate,  since  il  did  not 
h:i\-e  tia  larju.-  "irijwict"  spiki-  (;uul  suhseipieiit  "rinjtinjt")  measured  by  the 
fiiice  transducer  duriiiy,  the  short  period  when  the  water  column  is  impulsively 
ri'versod  in  tlie  iioiler.  Trecpieney  was  measured  with  a stop  watcli,  and  avera<;c 
hiiilcr  t eirpera  ture  hy  a t heri!iocou]>l  e recessed  into  the  boilei'  wall. 

TIk-  shutoff  \'al\e  location  was  such  that  the  stenm/water  interface  travelled 
Well  [ui.st  it.  So  haviiij;  taken  tfic  otlicr  measurements,  tlie  test  engineer 
closed  th.o  valve  while  the  interface  was  near  it?  maximum  travel,  and  then 
switclud  off  tl'O  Ik-. It  atul  w.-aited  for  the  cn,p,ine  to  cool.  The  water  tra])ped 
in  tlu-  duct  w.-is  siibs<.-ipienr  ly  dr.iincd  into  a metisuring  cylinder. 

Hecui-.e  the  addition  of  detergent  to  the  water  was  known  to  change  the  opera- 
tional ch.i  r.ici  er  i St  i e-;  of  p-ilsojets,  the  tests  wore  run  at  four  illfrcront 
concent  I, -It  i on  ratios  ( i nc  1 ud  i ii;’,  zero)  in  an  effort  to  see  whether  the  presence 
of  delorj’.ciit  modit'ieil  the-  ",51VR." 

The  "ct c.'dy  - St  at e"  boiler  temperature  variation  is  shown  in  Figure  22.  The 
data  border  on  statistic.il  r.'indomness,  but  a line  lias  been  drawn  because,  for 
many  enp.ines,  in  soni-  yeai-s  of  experience,  the  boiler  temperature  always  in- 
creases wh<in  tleterj'.cnt  is  .-idded. 

The  tlirin  t and  ju'riod  are  plotted  in  Figure  2,5,  :ind  are  seen  to  be  much  more 
regular  a.,  functions  of  concett  t rat  i on  ratio.  I!ut  tlie  stripped  water  volume, 
plott'-d  in  Figure  ."’1,  i •;  ajciin  very  scattered,  .so  that  the  suggested  reduction 
in  ".SlVI'"  with  i nc  r-.  .'i.  i lu;  concentration  verges  on  hypothesis. 

Not  counting  the  b.iilei-  cavity,  the  vo  I uiiie  of  the  Jiipe  "upstretim"  of  tlie  valve 
wa-.  2(1, in  ’ - .5-17. (i  ni.  liters.  Ihus  for  w.iter  alone 

”.SWK"  - 50/ .5  57,(1  = .089 

or  only  (lO"  of  the  estimate  based  on  ecpiation  (56). 


50 


BOILER  temper; 


550 


500 


DETERGENT  CONCENTRATION  IN  FARTS  PER  MILLION 


I'iK'n'i'  22.  Koilcr  teni]ieratiirc  as  a funct  ion  of  ilctc'rj^cnl  conccntrat  ion 
by  volume. 


THRUST  IN  POUNDS 


DETERGENT  CONCENTRATION  IN  PARTS  PER  MILLION 


I idlin'  27i.  Thrust  aiul  cyi.  Ic  I imo  for  the  iliict  iitiglcil  at  9“  (boilor  c'liJ 

lu'.'iM-  than  (1  i srha  rj’.o  oiul).  (Kims  ;t6()-2()7,  8/()/7.S).  -LSOb  watts 
Bell  mouth.  'TatiL  wat^'r  t cmiH'rat  uro  7!i°-V.S'’T. 


CYCLE  TIME  IN  SECONDS 


STRIPPED  water"  VOLUME  IN  ML. 


G'.'  WATER  AND  DETERGENT 

SALT  WATER  (3.63  PARTS  PER  MILLION) 


0 50  100  150  200 

DETERGENT  CONCENTRATION  PARTS  PER  MILLION 


^1.  Variation  of  ".st  r i j)|)i-il  watc'r"  voIuiik'  willi  ilot  I'l'gt'iil  concont  i .iLiun 
(Total  va  luo  i .s  ol)taini.Hl  by  aiMin,'.  i;il  to  fiituros  plot  toil  abovo, 
to  allow  for  fluid  whicli  ilid  not  tirain  out  during  iiu'asurc’nKnts . ) 


t" 


Tlu'ii.  arr,  of  i-,  •.■ra  1 f.iaiiM".  i-’lnali  n;  ly  aoi  mui  l fur  this  cl  i si' I'Cl’aiiay . 

(a.)  Tlu'  riiiicl  I'lii'.'.'  ' i liiini  1 s i v(  1 y reversed  by  the  lioiler, 
instead  el'  heiiiy,  s-teady  '-.tale.  So  the  exponent  n in 
Capiat  ion  (SI)  i.iay  have  a liiy.her  value,  especial  1\'  close 
t o t lie  hoi  1 er . 

(hi  I liiid  1 iirhnlenee  is  prohaMy  y.rcater  than  for  "nerinar' 
riciw,  hecause  of  that  indiicod  by  "jiiekiii;;  up"  the 
stripped  v.'ater  diirin;;  the  inr.troke  towards  the  boiler. 

'I'll  i s would  ahai  have  the  effect  of  increasing  n. 

(cl  .\  1 1 In Migh  iiovin:;  isoiv  slowly  than  the  interface,  the 

.strii'ped  water  is  j’.ciiei'ally  iiM-ing  outwards,  [larticu- 
larly  that  which  i close  to  the-  \'alve.  Thus,  the 
cpianlity  captured  should  he  soniewhat  less  ! ban  t he  v.  niu'c  t ed 
value.  This  effect,  coupled  with  not  c.U'sin.i  the  \'al\’c 
at  precisely  the  same  time  in  the  cycle,  nia\-  also  account 
for  some  of  the;  scatter  'n  figure  21. 

Having  said  that,  it  still  seems  lisely  (from  figure  .21)  that  the  presence  of 
detergent  further  increases  the  exponent  n to  give  a thinner  hoinidary  hn'cr. 

■fhe  jiresence  of  the  stripiu-cl  hoe'i 'n ry  layer  ele.srly  iiiodifies  the  derivation 
of  tiu;  eipiations  of  motion  from  c-iiwition  (1.2).  deferring  to  figure  .2.1,  W" 
assuiic  th.'it  tin  annulus  of  water,  of  thic'.kness  A,  is  left  behind,  and  Inteoiaes 
St. at  i unary . On  the  out.slr'oke,  at  an  interface  |to.silioii  x,  the  volume  of 
the  water  column  will  he 


OUTSTROKE  (dx/dt<0) 


INSTROKE  (dx/dt>0) 


lliiiil  foil  r I iMHM  t i i):i  . .•is'^iinicil  in  Ji-I  ormi  n i n;.’,  t lu-  oiiiiiitioii 
iii'ifion  with  .'1  -it  I' i I'j'cd  Ixni'i'l.iry  ];ivcf. 


Xj,  ■-  xll  - 2(VK)| 


I 1 - 2(^/U)  l|x\'  4 2(6/10  (X)-  I 


(60] 


p i.  i 11  still  liL-  re  1 .'It  c'd  t (1  ])  by  tHiii.i  t i on  ((.). 


I'or  tlvj  instrokc 


tiR 


(I\-  - l>) 


[2irR6x  + TtR  + iiR‘'xX| 


fl  - 2(6/R)](xx  - (X)- 


(C'l) 


ilK-  "blVR"  onrlicr  defined  !>>  otiiintioii  (bb)  )>ive5 
--  = "SIVR" 

so,  htisod  on  the  expor  itnent  ri  1 d;it;i 


2 


0.09 


1 - 2(6/R)  = 0.91 

so  that  while  metisurablc,  the  effect  of  the  stripiied  boundary  layer  is  not 
m;i  jor. 

The  complete  eqiiatiotis  of  motion  will  conttiin,  in  addition  to  those  above,  a 
tern  account  ini’  for  skin  fricti(jn,  which  will  be  discussed  later. 
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i 'U  iM  i :u)>!  i t .It  i Oil 

OlU'  V.-UI'l!  Il.lt  U’.lIlT  tl'  LOIll.lill 

I’ililv  I'i  ;;ii  . I',  v.iti.ii.',  iin-  >Mvfii 
I'or  ;i  Will’-  1 i'.i.po  1 .It  iin-  I'f 
I'Ik-  S'i  ! iih  i 1 i t y I'f  .\  i r i 
.N  i t rojp’ii 

Oxyjtcri 
(W-ii-tion  Itiuxiilc 


mil.  h di 
' ! ■ tvv  ; 

s solved  g,.'is. 

Some  tygiic.'il 

S-l'M' 

(.S^l- 

(1.  (bid 

(b  (bbl 

0.012 

(b  (Ido 

0.  (11  7 

(b  0106 

n.  (If).'! 

0.0. SI 

0. 018.6 

1 . S7 

0 . 9o 

0.2(> 

I'r.)!;!  till-  Ih'j;  i nn  i n;;  (if  [Hilsi-Jct  ilcvn.  )o|)!ncnl  , i iicoiiik-iis  i M e Rn;(  rclctii-c'  iliirinr, 
iHiiliiiP,  h.'i.s  |)(i.'".Cil  ;!  mr.iiher  of  pi-ot)  1 cm-. , nnj  the  (|Uii!it  i t i c;t  oJ'  y.at'  re  1 c.'i -icil  h.T\e 
'leemi  J ('.I'c.'itcr  rh.'iii  mij’.liT  be  siis|)cct(.(l  from  th.c  table  above,  assumipy,  tliat  oal\' 
air  is  (lissolvci!  in  tlie  v.’atcr.  Poo  1 -t'oi  1 i iii;  me.'isurcments  of  the  water  used 
ill  our  1 ahor.'it  ory  hav  shown  dissolved  "as  to  water  volumes  as  h i ;th  as  0.16; 

.'I  fij’.urc  not  i neon'-,  i St  (.'lit  with  the  bubble  volumes  sec-n  in  lipiirc  8, 

It  Seems  ele.'ir  that  dissolved  fttis  is  released  when  stenia  is  wade,  and  that  t!ie 
j;.-is  is  initi.'illy  mi  .red  with  the  steam.  Since  the  inteiface  nitiy  lie  very  irres'u- 
lar,  when  it  is  close  to  the  boiler,  it  seems  likely  tliat  at  least  some  of  the 
V.as  enters  the  water  at  this  time,  but  the  precise  mcchanisni  is  not  yot  known. 
()nce  r-omc  y,.'is  is  in  the  water,  in  the  form  of  a huhhlc,  it  oscillates  alone 
the  tube  with  tin-  surrounding,  water.  Dui-inp  this  oscillation,  faii'ls'  complex 
forces  act  on  the  huldile.  Its  size  ehanp.es  with  changes  in  the  local  static 
pressiifo,  ri'siiltiiig  in  cltaiij'.es  of  its  associated  virtual  water  mass,  ;ind  ;i 
viscous  drag  force  acts  on  it  as  well,  lixperiinenta]  ly,  it  seems  th;it  rlie  net 
effect  of  till  these  comple.x  forces  is  to  cause  the  bubble  to  migrate  towards 
the  exhaust  end  of  tlic  duct,  where  it  either  re-dissolves  or  is  expelled. 


rhe  e(|uations  governing  the  bubble  behavior  arc  of  some  interest.  Ilic  cC{uation 
for  the  radius  (rl  of  a stationary  biibhle  in  tin  inviscid  fluid  was  first  solved 
by  Rtivie  i gli . A niimher  of  workers  have  extended  this  solution  to  include 
surface  tension  tiiid  viscous  terms,  and  the  so-called  extended  Rayleigh  equation, 
in  ti  form  given  by  Scriven,^^  is 


who  11 ' 


r + d(dr/dt)'^  ( - ^ (dr/dt) 

, d or  or 


’’h 

i s 

ihc' 

huhhle  cavil)'  pressure 

>’1. 

i 

the. 

tiiihieiii  licpiid  pressun 

V 

i s 

Ml'.' 

I i (|ii  i d v i SCO  » i I y 

f) 

i s 

1 Ik’ 

siirf.ice  tension 

f> 


(6.1^ 


I 


It  X .'i’k!  :ir''  oriliiiatos  of  The  liquid  and  the  bubble,  its  motion  in  an 

in'.'iseiJ  fluid,  i .s  f.ovenu  d by'’'* 


d’xb 

' ‘ 

dt“ 


dr  ^'  "b 
r dt  dt 


(f’3) 


In  .1  viscous  fluid,  tliere  i-S  an  additional  term  due  to  bubble  drait,  v.hich  arises 
froin  the  movine,  bubble's  distortion  of  the  fluid.  If  we  use  the  formulation  of 
Moore'  ’ to  descrilie  this,  equation  fh.S)  becomes 


7 


.3  /dr 

V \dt 


+ 


l^vV’^b  ^ 

r |dr^  ^3 


(nr/2p)*^^ 


dx,  I 1/2 

-.tI 


54vi  dx 
"~T  dt 


(64) 


The  local  static  pressure  p^^  acting  on  the  bublile  can  be  derived  from  the 
cquation.s  of  motion  for  the  interface.  The  bubble  cavity  pressure  can  be 
calculated  by  assuming  isothermal  exp.-insion*  and  compression  so  that 


where 


Pb 


P„RT 


4 3 

trr 


UT 


K 

'"3 

r 


R is  the  universal  gas  constant 
T is  the  local  absolute  temperature 
K = .3m^RT/4rp  = 

ttij^  is  the  mass  of  gas  in  the  bubble 
is  its  radius  under  pressure 


(65) 


In  principle,  equations  (62),  (64)  and  (65)  enable  the  motion  of  a bubble  to  be 
calculated  if  the  motion  x = f(t)  is  known  for  the  liquid.  It  may  be  that  some 
pulscjct  configurations  could  have  an  oscillation  x = f(t)  which  would  cause 
bubbles  to  migrate  toward  the  boiler,  rather  than  away  from  it,  and  that  gas 
trapped  in  the  boiler  would  eventually  prevent  the  interface  from  entering 
and  steam  being  made. 

I-oi-  the  small  laboratory  ct\gines,  the  bubble  behavior  near  the  interface  ma.y  be 
dominated  by  the  turbulence  behind  the  interface,  an  effect  not  included  in  the 
foregoing  ccpiations.  Needless  to  say,  tlie  total  picture  of  bublile/water 
interaction  is  still  very  imperfectly  understood.  The  relatively  large  bubble 


* I’lc  discusses  the  conditions  appropriate  to  either  isothermal  or 

dialuitic  conditions  in  a bubble. 


r o 


voliriu'  ;i  1 so  "loans  that  the  water  slug  has  n fairly  low  sonic  velocity,  and 
should  pi'obaM)  he  wodelled  <as  a compress  i b 1 e twci-pliase  fluid,  at  least  during 
tlie  boiler  iaipaot  plia  .e.  High  speetl  movies  do  in  fact  show  that  a pronounced 
plane  shook  wave  travels  through  the  slug  after  each  boiler  impact. 

Although  in  principle,  the  equations  given  permit  the  motion  of  a bubble  in 
.an  incompressible  fluid  to  be  calculated,  they  ignore  another  effect,  which 
probabl)  has  an  important  influence  on  the  result.  This  is  the  finite  com- 
pressibility caused  by  the  presence  of  the  bubbles. 

In  place  of  the  simplicity  of  equation  (7),  inclu.sion  of  compressibility  terms 
leads  to  the  same  equations  as  those  used  in  the  study  of  "waterhammer" , ^ ^ which 
have  to  he  solved  numoTically  by  finite  difference  tecliniques;  woi'k  well  beyond 
the  scope  of  our  present  study.  Also,  as  may  lie  seen  from  IVi  jngaarden  ’ s^^  re- 
view, we  must  employ  a more  conijilex  representation  of  tlie  bubble  two-pluise  fluid 
than  the  ccimpress  i bl  c Newtonian  model  of  ivatcrhammer  analysis. 

One  manifestation  of  this  compressibility  can  ho  seen  in  I'ig.urc  2b,  which  sliows 
the  passage  of  the  pressure  wave  which  travels  along  the  duct  after  the  water 
intcrf.acc  has  impacted  the  boiler.  ,A  rarefaction  wave  follow.s  immediately 
after  the  sliock  front  has  passed,  in  which  the  bubbles  expand  greatly,  and 
others  seem  to  appear  from  nothing.  Possibly  cavitation? 

There  is  clearly  room  for  much  fruitful  research  in  this  area.  A possible 
way  to  start  is  to  separately  attack  the  calculation  of  bubble  migration  in  a 
representative,  externally  imposed  velocity  and  pressure  field,  and  to  solve 
tlie  compressible  fluid  colum'i  motion  problem  using  the  methods  developed  for 
wat orhammer  ana  1 ysi s . 

This  would  imply,  for  cxamiilc,  that  if  the  boiler  deceleration  length  were 
negligibly  short,  the  maximum  pressure  during  the  boiler  impact  phase  would 
be  limited  to* 

Ap  = pa(Xp  - (66) 

where  

a = 1/4^(1/K  + l)Cj/bb)  (67) 

where 

p = fluid  mass  density 

1)  = pipe  wall  thickness 

f.  Young's  inudulus  for  tlie  pipe  material 
I)  = pip*-'  diameter 

K = hulk  modulus  of  the  fluid 


* See  pp.  30  of  Peference  29. 


r-'igure  26.  These  sequenced  frames 
from  a high  speed  movie  show  a 
low  pressure  wave  moving  from 
left  to  right,  as  exemplified 
by  the  enlarged  bubbles  or 
cavities.  Film  speed  was 
approximately  1000  frames/scc. 


Cj  is  a constant  u-hicli  cI(.']>LMids  on  details  of  the  I)ilie  fixation 
■ 5/1  - a for  a piilseiet 
= Poisson's  ratio  for  the  pipe  material 

3 

Kt)uj;h  order  of  maj'.nitnde  figures  are  p = 1.8  sliigs/ft‘,  x.,  = 20  ft/sec, 

= -4(1  ft /sec,  a - oOO  ft /sec. 

A|i  - .■=;2.400  IS.-'S  p^ 

wliich  is  certainly  adecpiate  for  good  t hermotl>naii:i  c efficitncy.* 

There  is  clearly  a need  for  more  detailed  ex|)eri menta  1 observ'at ions  of  bubble 
behavior.  Some  can  always  he  seen  emerging  from  the  duct  exit  during  the  out- 
stroke;  and  subjectively,  the  shortt'r  the  engitie,  tlie  greater  tlie  volume  of 
hubbies.  On  the  other  h.ind.  tlte  reliability  of  engine  operation  can  be  improved 
by  placing  "bleeders"  in  the  boiler  attachment  flange.  These  are  sm.al  1 j’ij'cs, 

1/8  incli  1.1).,  which  m.iy  either  project  'tito  the  engine  duct,  or  be  flusli  with 
tlie  wall.  The  optimum  location  for  these  bleeders  will  be  discussed  later,  but 
in  this  condition,  they  di scharge  gas  bubbles  and  engine  ])erformance  is  improved. 
This  would  seem  to  indicate  that  some  gas  does  not  naturally  migrate  away  from 
the  boiler,  but  stays  close  to  it,  and  somehow  impedes  boiler  ])enctr;it  i on  by 
the  water. 


Th I . H’ e c^  of  lie t erg e n t 

At  an  early  stage  in  pulsejet  development,  we  were  faced  with  a momentarily 
baffling  jirohlem.  A pulsejet  liad  been  running  satisfactorily  for  months,  but 
a newly  built  one,  of  identical  design,  refused  to  function  at  all.  It  was 
eventually  decided  that  the  internal  surfaces  of  the  new  one  must  be  greasy  or 
in  some  other  way  contaminated  by'  the  various  forming  operations,  and  that  this 
was  the  reason  for  the  difference.  (Rhodes  and  Bridges^'^  found  that  the  rate 
of  boiling  of  water  on  steel  was  considerably  reducctl  by  a trace  of  mineral 
oil.  They  also  found  that  a little  sodium  carbonate  could  change  film  boiling 
to  nucleate  boiling.)  Pouring  detergent  down  the  engine  enabled  it  to  start 
running  weakly,  and  gradually,  over  a period  of  an  hour  or  so,  its  performance 
came  to  match  that  of  its  fellow. 

It  was  natural,  then,  to  see  what  effect  detergent  had  on  the  obscrx’ahle  opera- 
ting parameters.  The  steam/watcr  interface,  for  example,  was  no  longer  sharply 
defined,  hut  had  a frotliy  "head"  on  it  and,  as  we  have  already  seen,  the  addi- 
tion of  small  cpiantities  of  detergent  seems  to  reduce  the  thickness  of  the 
"stri|iped  boundary  layer,"  and  causes  some-  variation  in  thrust  and  frequency. 

In  I'igurcs  .'7  - 29,  we  present  similar  data  for  a much  shorter  engine  which 
runs  at  tihout  .1  Hertz.  The  variation  of  thrust  with  concentration  ratio  is 
seen  to  he  quite  complicated,  and  deitends  on  the  angle  of  the  duct  To  the 
horizontal.  This  angle  ;ilso  influences  the  boiler  temperature,  but  not  the 
operating  frequency.  As  indicated  in  Tigure  28,  there  is  initially  a marked 

* lor  .'in  air  breathing  pulsejet,  assuming  sonic  velocity  in  each  direction,  the 
same  analysis  gives  Ap  ' 2yP<..  " 2.6  p„,.  This  low  value  explains  the  poor 
efficiency  of  air  breathing  piilsejets. 


BOILER  TEMPERATURE  ®F  THRUST  IN  POUNDS 
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DETERGENT  CONCENTRATION,  PARTS  PER  MILLION 

liHiirc  27.  The  cffci.-t  of  dctcrsenl;  (Tide)  on  thrust  :iiul  boiler  t enii)erattire , 
(I’PM  by  weiy,)it).  Duct  angles  of  - 11 " and  -1°  mult  i -hole  boiler. 
2.7.S  ft  duct  length. 
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STROKE  IN  FEET  CYCLE  TIME  IN  SECONDS 
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ri'.’.uTe  ^8.  The  efFeet  of  iletcrgeiit  on  eyelo  t ime  ;inil  sTrt'ke,  for  the 
conditions  of  l■ij;llrt•  71. 

6.1 


DYNAMIC  PRESSURE  -g- /}V  IN  LB/IN 


CM 


I'iguro  29. 


v'^iipori niposcd  dynamic  pressure- time  histories  from  the  experi- 
ments of  Figures  27  and  28.  (Zero  time  for  each  case  is 
offset  0.1  sees  for  clarity).  Runs  243-248  (7/31/7.3). 
Multi-hole  boiler;  Temperature  275-310°F.  Duct  Angle  -I®. 
Duct  I.ongth  2.75  ft. 


(Note:  The  change  in  slope  passing,  through  zero  is  thought 
to  be  due  to  tlie  finite  frequency  response  of  the 
transducers.  Differenl  transducers  were  used  for 
instroke  and  outstroke,  since  a pilot  head  is  uni- 
d i reel i ona  I . ) 
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incrc.isc  in  stroke  (X  - to  about  50  PPM.  Stroke  is  then  constant  to 

90  PP:'I,  after  which  it  decreases  so  markedly  that  it  does  not  emerge  into 

the  ti  ansporent  poi'tlon  of  tlie  duct.  The  dynamic  pressure-time  histories 
plotted  in  Figure  29  also  show  significant  variations. 


No  data  is  sliown  for  clear  water  because  this  particular  engine  would  not  run 
continuously  without  some  detergent  being  present.  "Full  power"  (6156  watts) 
was  insufficient  to  liold  the  boiler  at  working  temperature  using  clear  water. 


We  offer  these  observations  without  attempting  to  explain  them.  The  detergent 
may  affect  tlie  heat  transfer  coefficient  in  both  boiling  and  condensing  phases. 
It  may  also  improve  the  water  column's  "capture"  of  gases,  so  that  migration 
of  bubbles  away  from  the  boiler  is  more  reliable.  The  fact  that  tlie  boiler 
rccjuires  less  heat  injnit,  for  a given  temperature  when  detergent  is  present  is 
possibly  connected  with  the  frothy  "head"  at  the  interface.  That  is  to  say, 
loss  solid  water  penetrates  the  boiler,  so  that  less  heat  is  carried  away 
from  its  walls.  Alternatively,  or  perhaps  concurrently,  the  presence  of  the 
detergent  in  the  water  may- reduce  the  heat  transfer  coefficient. 

In  experimental  boiling  of  n-pentanc,  Berenson^^  found  that  a clean  oxidized 
surface  had  the  same  pool-boiling  characteristics  as  an  unoxidized  one,  but 
that  exposure  to  the  air  greatly,  increased  heat  transfer.^  Me  also  found  that 
the  addition  of  a surface-active  agent,  such  as  oleic  acid,  greatly  increased 
heat  transfer.  Thtis,  there  is  precedent  for  our  findings. 

nerenson  believes  that  oleic  acid  (or  a dusty  surface)  influences  pool  boiling 
because  tlie  contact  angle  is  changed  from  10°  (for  n-pentane  on  clean  copper) 
to  spreading;  i.e.  no  contact  angle. 

A Comparison  Between  Theory  and  Experiment 

The  impulsive  ])ulsejet  cycle  equations  presented  earlier  are  compared  with  an 
actual  velocity  time  history  (from  pitot-static  measurements)  in  Figure  30. 

This  comparison  was  made  by  using  the  observed  stroke  and  frequency  to  deter- 
mine the  equivalent  constant  pressure  P from  equation  (41),  and  aligning  the 
theoretical  and  experimental  curves  at  the  time  (t)  for  x = x . Because  skin 
friction  is  neglected  in  the  model,  the  agreement  is  worst  on*"  the  outstroke, 
when  the  velocity  is  highest. 

The  general  equations  of  motion,  including  a conventional  D'Arcy  friction 
factor  (f)  term,  liave  been  programmed  for  the  digital  computer,  and  in  Figure 
31,  the  output  of  this  program  is  compared  with  the  same  experimental  data; 
again  assuming  constant  pressure.:  As  for  the  simpler  analytical  model,  it's 
assumed  that  the  boiler  is  a "black  box"  which  accept.s  the  interface  entry  at 
velocity  Xj,  and  instantly  expels  it  at  x . The  latter  velocity  would  usually 
be  computed  from  an  assumed  or  calculated  "boiler  efficiency"  (defined  by, 
equation  51)  but  in  this  case  x^^  Was  chosen  as  the  value  which' gave  the  correct 
stroke  and  frequency,  in  conjunction  with  a constant  pressure.  . Agreement  with 
experiment  is  seen  in  Figure  .31  to  be  much  better;  possibly  as  good  ns  can  be 
expected  with  the  constant  steam  pressure  assumption,  and  at  this  stage  in  the 
evolution  of  the  technology,  ; ; - : , i v-  : -v 
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includes  skin  friction  ter" 


► 

i'.'o  K'luiic  t lu'  i iKi>iii|)ri-':-^  i h 1 1-  t'liiiil  i iii’ i .'i  1 1 )'  s.'i  t i ■.  I'.i  r t (u  y , 

' .lurin;,;  t hr  hoi  Irr  iiii)).ii  1 rvrnt  , ulu'ii  il  t'r  i,(  ci  ■ s.n  y In  ! rr;jt  if 

:r  , I' ifss  i h I . lloKinri-.  i l ' ilrnhfriil  i.  Ini  ir' r tin-  i .m,'.  ■.  ;it  i oii-i  1 Ii'\ri\  '-kin 
frirlMiii  rr  1 .1 1 i rush  i |i  i •,  \rr\’  n- 1 rvaii  1 lo  tin'  r'-i  i ! I i i ay  flow  if  a 

inilsrirt.  (Sa-aI^'  shr.'-.ii  ! f,  i ■'  I'or  I aia  i na  r'  I'lnv,  '.•tii.'.  ! hr  flaiil  i rr-r  i ] ) a i i-l 

i -o  ' ,ia  1 1 y . ) 'lln-rr  is  llui'c-ri>rr  a ar-r  Hu'  min;;  I h.r  ■.  j ii| ) ] ]•  la- 1 a i i on  sli ; |i 


' *'ri'  I rt  ion 


K|,  ,.(x)' 


\>liero  K..  is  (.Ictr  niri  nod  from  oxpri-ii'rnt  ; perhaps  with  an  apj'uratus  based 
"drinkini;  straw”  prol)lcm  solved  'oarlier.  Also,  win  a the  diic;t  Iriijjth  is 
oo.iiparctl  with  tltc  stroke  (as  is  tlio  case  with  water  juilsejct  pomps)  x nia 
replaced  by  a constant  x (the  mean  value)  whiih  then  permits  the  closed 
"inviscid”  equations  to  1)0  used  to  obtain  solutions. 


Ion;; 
y be 
form 


(>8 


Tllli  rill.kMODYWMICS  Ol-  W)(I.|;K  IMPAC'I' 


Pj  ilcr  lU:;' i TiiAc 


In  nrclc-r  to  clot  ot-’,i  no  the  he.it  t ran-,  rcri-ed  to  tlic  hoilo)-,  ue  liave  to  cstalilish 
a physio, 'll  r.uulc!  of  tlu'  l>oi  lor  i)enet  rat  i on , and  of  coui’So,  v.o  have  >'C't  to 
ohsorve.'  tills  in  the  1 .'iho/'a tory . 'I'he  simjilest  possible  model  assumes  tliat  the 
interface  will  renain  i)lanaT  in  the  boiler,  and  that  its  motion  is  arrested 
by  compiess i on  of  rtic  steam  and/or  {'.as  trapped  above  it.  But  this,  comi'rcssion 
process  is  ver)’  complc.x,  since  the  gas/vapor  is  initially  heated  by  the  liot 
walls,  and  r,;a>',  at  high  [iressure,  reject  heat  back.  So  at  the  pre.sent  stage 
there  is  perhaps  little  to  be  gained  by  using,  a pressure  re  la  t i ons'n  i ji  more 
comple.':  tb.an  Ayi  - K ( .x  - ,\n). 


.Since  tbe  deceleration  lakes  place  within  tlie  boiler,  the  x’aliic  of  .x  in  equation 
(7)  is  not  goiiK;  to  vary  much,  and  can  Jiist  as  well  be  written  as  .a  constant, 

X|^.  he  also  assume'  that  since  (xx)  is  large,  the  velocilx'  squared  term  in 
equation  (7)  can  be  neglected.  (This  iiii[)lies  that  no  noc.le  is  fitted.)  So 
tlie  equation  of  motion  in  tlio  boiler  simplifies  to 


or,  writing  z = x - x^,  z 


etc. 


• dz 
dz 


(6S) 


1 ,'2  • 2-  Kz 

2 - "l5  ^ = 2^x-„ 


(69) 


SO 


= + Kz^/pXj^ 


N'ot  e 


that  7. 


0 when  z 


z , the  maximum  value,  so 
c 


K 


-PXi, 


Si ) 


(70) 


69 


56 


and  Ti-  s idcnci.'  t i 


t 


R 


sin 


f71) 


f72) 


If  the  heat  tranaf'tr  coof f icit-nt  and  the  teia])eraTnre  difference  AT  (wall  to 
water!  are  constant,  the  heat  transferred  in  one  penetration  is 


2'!fRhAT(t^/Z|^ 


sin  0 do  = 


for  ,j  typical  engine  - 20  ft/sec,  = (i.4  ft,  AT  - 
and  t lie  exiievimontal  ly  d'  termined  heat  transfer  is  1.65 
this  case 


lirRh.ATfzJ/^p)  (73) 

220”r,*  R = .0286  ft 
Btii/piil  so.  For 


t 


R 


0.4  II 
■"2(i" 


.063  secs 


1 . 65  X 20 

' " dVV.o^SO  ■x'220  7a  TnT 

(0.4)“^ 

= 2.61  Btu/ft^sec.'’f 

- 0391  Btu/fr^hr.”!- 

Actiril  ly,  we  .shall  sngi;e,sl  later  that  the  water  is  only  in  contact  with  the 
boiler  w.all  diirinj;  the  instrokc,  except  iiossibly  for  that  near  the  boiler 
entr.incc  when  the  interface  is  near  maxiinuni  penetration.  Also,  an  insulating 
steam  layer  may  be  present  for  tiic  earlier  portions  of  the  instruke.  If  these 
hypotheses  arc  coirect,  then  the  true  value  of  h may  be  significantly 
higher.  Finally,  the  assumed  boiler  penetration  of  0.4  ft  may  be  too  high. 

If  it  proves  to  be  less,  this  would  further  increase  the  required  value  of  h. 

2 

.Nominal  values  of  h up  to  12,000  Btu/ft  hr.'’F  have  been  observed  on  some 
engines,  ainij  this  is  considerably  greater  than  is  usual  in  the  heat  transfer 
literature  for  pool  boiling.  It  therefore  is  important  to  understand  why 
this  should  he  so. 


Heat  Transfer  fsl imates 

Ke  start  with  the  assumiition  (to  be  discussed  more  fully  later)  that  water 
is  in  direct  contact  witli  the  boiler  wall  duriitg.  penetration,  ihen  a conven- 
ient way  of  computing  li  is  to  use  the  Reynolds  analog)'.  This  is  in  agreement 


* The  lioilcr  ttmip'  rature  was  .STO^F.  We  Itore  assume  150'’F  for  the  entering 
water,  but  it  could  be  anywhere  between  13t)''l-  and  20(i°F. 
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ii  is  r lu-  II  i|u  i il  V i siM  ' i t > 

n is  ; III-  snrC.K't.'  ti.-iisi<Ki 


v\' i t il  ' ■ \;h' r i I ."i  t I '.iita  for  fluids  ( sur!)  us  ivator)  li.'ivin;'.  I’raiiiltl  mini'ieri 
un  i t V , ' i . o . 

p - 1-0 

r K 


kIic  rc 


is  rhi.‘  specific  lic:if  at  constant  prcssin’c 
ij  is  tlic  ahsoliitc  (or  dynamic)  viscosity 
k is  tlie  thermal  conductivity, 
iiie  Keynolds  ;malop,y  leads  to 

St.  “ * (the  Stanton  number) 

pCpU  8 


Kliere 


f is  the  conventional  D’Arcy  friction  factor 
=•  Ap/i'.u-(L/l)) 


p is  the  mass  density 


u 

It 

L 

It 

Ap 

For  a smooth  duct, 


is  the  averape  fluid  velocity 

is  the  hoiit  transfer  coefficient 

is  the  I’ipe  lenpth 

is  ttie  pipe  diameter 

is  the  pressure  lo.ss  over  tlsc  Icnptli  1. 

tiic  Nlkuradse  ec|uatlon  gives 


1, 

/f 


2 logjjjCv'f  Re)  - 0.8 


S i nee 


Rc 


uO  _ iiDp 

V 1) 


III) 

'Hie  variation  of  f 


"fp 


(Re  f) 


Re  with  Rc  is  plotted  in  Figure  ,^2. 
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tu-ar 


(7-n 


I 


*■  Sec  pp.  3n  of  foroiico  29. 


r.9 


OrluM-  oiiiinf  i oils  ill  u i (i'.'Sjiro.’id  use  foi-  l'i)i\eil  convention  turlnilent  pi))',  floiv 
are  i ho  foK  urn  eauarion,  the  Sicdc-r  Tati-  ciiuat  i on  , aiul  the  1)  i 1 1 us- Boo  1 1 or 
rel  at  Loiish  i :'■* 

= n.  02.3 (Re ’ 

Siihst  j t lit  i no  k - 

liO  - 0.023  C|,)i(Ro)"-^fl'^)  (7.3) 

Th  is  is  also  plf'tteil  in  figure  32,  as  a check  on  the  analy.sis.  Ai’reein'. ait  is 

scon  to  be  verv  i;ooil  when  P = 1.0. 

' ' r 

.Now  for 

Temperature  - 60°I-'  200'’r  300“f  lOO^r  410°r 

iiC|,  (Rtii/hr.  ft . °l' 1 ■-  2,7189  0.7418  0.4()23  0.3.396  ().33(i2 

At  the  temperatures  associated  with  a .SWIM  boiler,  uCp  v 0.4  therefore.  Then 
for 


Re  = 4 X 10 

lO'’ 

10''’  10^’ 

10^ 

hO  - 7.8 

13.0 

88 . 3 580 

4030 

(Btu/hr.ft^°F) 

Now  consider  201)'’ 

H water  entering  a boiler 

;it  20  ft/soc 

For 

1)  - 

0.  1 

1.0 

10 

inches 

Re  = 

4.8  X lo'’ 

4.3  X 10^ 

4.8  x lO*-' 

f Re**  = 

9.7  X 10^ 

6.2  X 10'^ 

4. 1 X 10^ 

hi)  - 

48.5 

310 

2050 

h = 

3820 

3720 

2460 

Btii/hr.  ft 

* The  suffix  "b" 

refers  to 

the  bulk  of  the 

fluid. 

**  l'akin;>  f for  a hydrudyrvnnieal  ly  smooth  wall.  Practical  boiler  walls  are 
roup.her  than  this. 
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'IIks..'  boat  aro  of  tlio  samo  oialor  as  Vvc  oliscrvo  cxporinoiital  ly. 

C’rrati'i'  t lian  nor-'.al  t urtni  I once , duo  lo  tho  strippod  boundary  layer  ivater,  inny 
oxplain  wli_\'  til'.'  oxiu- r i i;;0ii  t a 1 values  aro  sciRii-v.-hat  liijdior.  Note  also  that  f 
inoroasos  uith  sin'fao<-  taaij-hacss ; by  as  miioh  a;',  an  ojalor  of  i:ia»nitudo  at  lary.c 
Uo.  ^ ' Thus  tho  boat  tiauisfor  in  a f.ivp.l  hoi  lor  is  aclocpiatcly  oxplainoJ  by 
exist  inr;  t!;oo?'.v,  onoi-  i.c  rocojpiixo  that  stoaiii  uill  not  form  until  tl'o  wator 
coliriin  !;  IS  roboin  led , booauso  local  prossuro  is  liiaher  than  tlio  vapor  prossiiia'. 

lloKovor,  ti  cautionary  noto  should  ho  sounded.  In  a typictil  ono-iiuli  diameter 
])ul.seje-,  th'  lieat  transferred  per  lioiler  iRi|)ai:t  i .s  of  the  order  t'f  2-.T  !Uu. 

If  all  tills  i.are  used  to  make  stenni,  the  steam  wei  ght/piil  sc  would  he 


Q 

11 


2 . r. 
lyob 


.1)021  Ih. 


IVh  i 1 1'  in  the  ei'cipressed  1 tcpiid  sfap.i',  the  cot  res()ond  i np.  rolumc  of  water  is 
O.of,  in-"’.  Tliis  is  small  corijiared  with  flu'  total  volume  of,  say,  I in-’’  in  u 
1"  \ straij’lit  boilei  . In  f.aet  , it  r(.-pr<,'S<,';its  an  outer  atinulus  v.hieb  Ls 
oiilv  four  thoie-aiidt hs  of  an  inch  tliick.  Due  to  losses  in  heatin;',  water  wliicli 
is  not  turned  to  .steam,  the  annulus  of  "useful”  water  may  l<c  even  thinner 
than  this;  porhtips  less  tliin  a thousandth  of  an  inch.  Thus,  conventional 
(cjuas  i - stilt  i c ) ide.ar.  of  heat  tnui.sfer  may  sometimes  mislead  us. 


Tl^^'  l!oi  1 er  Cyc  1 e 

As  a reference  point,  it  is  instructive  to  first  ct'iisider  tho  .Savory  puiiiiiiii;.; 
cycle,  in  which  writer  is  e.xpclled  from  a vessel  hy  steam  p-os.sure,  tin.i  a fresh 
supjily  drawn  in  tis  the  steam  condenses.  The  idealized  cycle  for  tins  is 
presented  in  Tiipire  S.i,  and  tlie  cffieiency  e<a!>  he  obtained  directl)'  from  the 
ste.im  t.ihles.^'^ 


hot  II 

V 

.S 

A 

fd 

.1 

The  amount 
1)  toe  in 


--  steam  enthalpy,  lUu/llt. 

steam  specific  volume,  ff/lh 

- total  stroke  or  travel  of  the  water 

- duct  cross-sectional  area 

= steam  pressure  during  disehtirge 

steam  pressure  during  induction 

mecli'inieal  e(|ui  valent  of  heat 

of  steam  rccpiired  to  discharge  the  water 
Pigure  3.'i)  is 


through 


the 


stroke  ? 


(from 


1V„  - AS/v  Ih 

S o 
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111  - (.\S/v  )ll  - II'  ) 

kIi.ti.’  II‘  is  till'  intli;ilpy  of  v.-.-ifov  nt  an  apjii  opriato  1 eirij»o  ra  ^u  ro ; oi  flier  that  of 
the  eoiidensiiii;  steaia.  or  that  of  --oiik'  oatsiile  supply  source.  In  pr.Klice, 
tluM'c  is  little  (I  i ffcr.'iK  e between  tlu-  two,  aiul  hijMire  is  drawn  for  the 
first  ease. 

Since  the  u.scftil  work  done  per  cycle  is 

K = AS(p^j  - i.j) 

the  efficiency  of  the  cycle  is 

h.  ■ ''l  ' 

'b  ■ JiL  ■■  jCIT'-iT*! 

1 o 


,As  shown  in  I'iyure  .^4,  the  efficiency  of  this  cycle  is  quite  low. 

Now  consider  the  I’ankine  cycle  in  I'iyiire  S5,  whore  steam  is  iiiatle  or  admitted 
for  only  jinrt  of  tlu'  cycle,  and  expands  to  do  tidditional  work.  (Tins  cycle 
was  first  employed  hy  .James  Watt).  For  the  water  to  be  completely  discharged 
without  kinetic  energy  in  the  discharge 

I (p  - tiv  = 0 AS<f(p^j  - P^)  + A / (p  - pj  dy  (7b) 
%iiN 

As'^iining  isothermal  cxjiansion 


P = Po(<l<S/y) 

Making  this  substitution  in  equation  (76) 


'{>‘i(P{,  - P„)  - (1  - <t’)Sp^  ' p 


l/<!’Il  + log  {1/-HI 


(P,./P.,.)/|1  ^ log  (l/'Ml 


7(, 


UC.IC.K(]tNI  tUNUfcNIKAI  lUN,  HANTS  PtK  MILLION 


f’i'TUiC!  :'8.  Ttu'  (.'fFcet  of  dftcrefut  on  lyclc  tinu'  :ini.l  srrokf,  for  the 
c'C'mlitions  of  I'ijiore  ?7. 

f.o 


SUCTION  PRESSURE,  LB/ IN*  ABS 


(•ij.’,ure  81.  Idi’al  Savei’y  ffricieney  for  Ket  Sieaiii. 


6 


1 l;Mlic 


I cloa  1 i :■  <il  Statir  RaiiKino  l’iii;i|)  i ii;; 


I '/I  J . . , 

cxfic  r i iii'Mif  is  sec'ii  in  SI  to  In*  mucli  ticttcr;  possibly  :is  good  as  can  be 

expecti'd  witli  the  constant  steam  pressure  assumption,  and  at  this  star.o  in  the 
evolution  oT  the  technology. 
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llu'  useiu! 


V +.\s 
MIN'  MIN 

J !'  j n>  - i>„ ) + j (IV 

CVCI.f^  ^VllN 


- P.ldV  - ASfp  - p ) 


V + AS 
MIN' 


(as  before) 

(80) 


Ihe  ai.ouiit  of  steam  re(|uircd  is 


V. 

S 


A(}.S 

V 

o 


lb 


The  heat  re<]uired  is  Ihcrefore 


A!)|S 


(II 


II*) 


So  the  efficiency  is 


jaTsoT  - iFj 

o 


V ( p • p . ) 
O ‘Jj 

'<))J  TH^  - II*)- 


where  II  is  now  based  on 
o 



’’o  " <{.[r“*  log  (i/T)i 

(y  - 1)  IV 

'(’(Y  - '^) 


for  isothermal  steam  expansion 


for  adiabatic  steam  expansion 


where 


Y 


the  ratio  of  specific  heats 
1.1  - 1.2  for  wet  ste.am 


(81) 


(82) 


As  shovai  in  I'igure  7i0 , this  permits  siibstant  i .i  1 ly  higher  efficiencies,  although 
;it  the  cost  of  higher  initial  ste.am  temperatures. 

The  cycle  efficiency  of  present  day  water  pulsejets  is  substantially  lowo’-  than 
these  values,  p resiimab 1 y due  to  the  heat  lost  by  the  steam  to  the  duct  walls 
.and  the  stri|)pfd  boundary  lave)-. 
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THERMAL  EFFICIENCY 


INITIAL  STEAM  TEMPERATURE  IN 


When  the  v.'ator  column  contacts  the  boiler  wall,  it  acquires  heat  from  the  wall 
at  a very  much  higher  rate  than  the  external  heat  supply.  Thus,  the  wall 
temperature  falls  as  that  of  the  water  rises,  as  indicated  in  Figure  37. 

Tn  effect,  the  thick  boiler  wall  acquires  and  stores  heat  while  the  liquid 
interface  is  elsewhere,  and  that  portion  close  to  the  surface  discharges  heat 
"capacitively”  to  tlie  water  when  it  renews  its  contact  with  the  boiler  wall. 

The  Iieat  acquired  by  the  water  is  not  distributed  evenly  across  the  column 
(unless  its  turbulence  is  very  large)  but  is  mostly  retained  by  the  water  close 
to  the  boiler  wall,  as  indicated  in  Figure  37.  (Note  that  the  water  need  not 
vaporize  during  the  instroke  but  water  at  a given  temperature  will  flash  to 
steam  on  the  outstroke.  when  the  pressure  has  fallen  to  the  appropriate  value.) 
Since  (as  we  shall  prove  later)  the  distance  over  which  the  temperature  varies 
appreciably  (x'  in  Figure  37)  is  small  in  relation  to  either  the  radius  of 
the  boiler  or  its  wall  thickness,  we  can  use  as  a one-dimensional  model  two 
semi- infinite  slabs.  Unfortunately,  even  for  this  simplified  model,  wo  cannot 
precisely  calculate  the  temperature  profiles  without  recourse  to  finite 
clement  numerical  analysis,  oecause  the  fluid  is  moving  in-the  y direction, 
relative  to  the  boiler  wall. 


At  each  boiler  wall  location  y,  we  can  estimate  the  wall  temperature  from  the 
relationship^^ 

2 

T = T^  - (Tj,-T^)[erfc(n)  - erfc  (n+B)] 


(8.3) 


where 


erfc(q)  = J o 

'q 


CO 

cix 

J n 


1 - ^ 


i 


dX 


k/Cppg 


n = x/2»/aJ 

3 = hv^/k 

a 

= initial  wall  temperature 

bulk  temperature  of  the  fluid 

wall  temperature  at  some  distance  x from  the  surface 
film  heat  transfer  coefficient 
elapsed  time 

conductivity  of  the  wall 


'F 

T 

h 

0 

k 


% 1/ 


03¥  AVAILABLE  COPY 


CONSTANT 

TEMPERATUI^E 

EEFOi?E 

A:lF?iVAL 

OF  THE 

V/ATER 

COLUMN 

JUST  AFTER 
ARRIVAL  Or  TH 
INTERFACE 


AFTER  THE 
INTERFACE  IS 
WELL  PAST 


\ .1  r i ,1 1 i III!  >>l  lioilcf  iv;il!  iiiul  ivati-r  t I'niK' r,i  Mi 


specifii.'  boat  i)f  t(u'  wall 


,o  - s]>t.c  i f i 0 Weight  of  tlu'  Vv-all 
riu'  correspoiul i ng  stress  in  the  material  is  approximately 

<■  ■ (T  “'si 

wlicre  f.  = Young's  modulus  for  the  material 

e -■  Poisson's  ratio  for  the  inal.i>rial 

o - the  theiT.ial  expansion  coefficient 

Tfic  greatest  t C'lnporature  cliange  and  therefore  the  greatest  stress  occurs  on 
the  surface,  where,  from  ecpiation  (<S^) 


T - '1' 
o S 


orfc  (;^ ) 


(SS) 


This  function  Is  jilotted  in  l•.igtlre  38.  If,  for  example,  wc  lake  the  numl)e;^\s 
used  earlier  to  determine  the  heat  transfer  coefficient  of  h = 57.36  Btu/ft"hr. ‘’F, 
and  a total  residence  time  of  .063  seconds,  wc  find  that  at  .0313  seconds  fwhen 
the  inward  motion  is  stopped) 

6 - 0.103,  T - Tc  - 220  x 0.103  = 22.7  '’F 

o S 

f 4762  Ih/in^ 

max 


a stress  figure  v.hich  is  not  intolerable.  From  I’igurc  38,  factors  tending  to 
increase  tlie  thermal  shock  arc 


higlicr  heat  transfer  coefficients 
loiigei'  rosid(Miee  t ime 
.smaller  values  of  /pg  (!|,k 


(\o’-e  that  for  Copper  »'pgC|,l 

.pno.-i 

1 / 1 a 

hr*'  • ft  “F/Ptii 

Hr.iss  " 

.('IS 

,,,.l/-’,-t^'’‘'F/lUii 

A Ill'll  ilium  " 

.017 

,,,,1/dp,  2oi,/|^tu 

■StaiiilesM  .Steel  " 

~ ,01' 

,,,,l/2pt 

... 

iS3 

0 0.5  LO  1.5  2.0 

SS.  Rnilfi-  surf. in-  t '•liiptratiiro  ) .'is  a fnaclion  of  t ho  )in  raiiiotc- r 
V , froi'i  I'uu.if  i oi’  ( Sh  ) . 


81 


31' Tenpc r;.i t lu' c Distribution  in  the  Boiler 

Since  the  change  in  wall  tempciaiturc  is  only  22.1°T-  (out  of  a total  AT  of 
22S‘’F)  a first  approximation  to  the  temperature  distribution  in  the  water 
is  obtained  by  assuming  the  boiler  wall  temperature  to  be  constant.  We  can 
then  employ  equation  (83)  if  wc  know  the  thermal  diffusivity  (a)  of  the 
water.  For  laminar  flow  of  hot  water  (Tp  = 200°F) 

k ~ 0.39  Btu/ft^hr. °F/hr 

Cp  = 1.0  Btu/lb 

pg  60.0  lb 

a ==  .0063  ft^/hr 

In  this  model  of  the  water  heating,  the  highest  temperature  occurs  where  the - 

interface  water  is  in  contact  with  the  wall.  Here  the  pressure  (from  equations 
68  and  70)  is 

Ap  = pXp  (zp/z^)^z  (86) 

where 

z = X - Xp,  the  boiler  penetration  at  time  t 

z = the  maximum  boiler  penetration 

0 

» 

Zp  = tlie  velocity  at  which  the  intoa'facc  enters  the  boiler 
From  equation  (71) 


where  the  suffix  IF  denotes  the  interface. 

At  the  point  where  the  interface  has  been  arrested,  we  have  from  equation  (72) 


2 ("c''’"!)’ 


A point  in  the  column  which  is  a distance  q behind  the  interface  has  therefore 
been  in  contact  with  the  boiler  wall  for  a duration  of 


t 


q 


(z^/Zp)  [(Tf/2)  - sin"^(q/Zj.)] 


(88) 


BEST  AVAILABLE  COPY 


Usinr,  tlris  relationship  with  equation  (Sa)  , we  can 
clistvi  hut  ion  tliroughout  the  water  at  the  Instant  of 
matter,  during  the  subsequent  outstroke. 


calculate  the  temperature 
arrestment;  and  for  that 


In  I'igurc  39,  for  the  pulsejct  example  previously  used  in  this  section,  we 
havi-c  plotted  tlic  water  surface  temperature  at  the  interface,  using  equation  (87) 
to  determine  time,  and  equation  (85)  for  the  temperature.  The  corresponding 
critical  pressure  is  read  from  the  steam  tables.  The  simple  linear  pressure 
relationship  for  actual  pressure  (on  which  equation  87  is  based)  is  also 
simwn.  Insofar  as  this  is  correct,  the  model  indicates  that  some  boiling  will 
occur  just  after  entry  because  > P-  Then  heat  transfer  will  be  reduced, 

lierhaps  to  a negligible  value  and  will  commence  again  when  ponjo,  < P> 

and  water  is  again  in  contact  with  the  boiler  surface.  In  Figure 
39,  wc  have  not  attempted  to  model  this,  hut  have  assumed  complete  water  con- 
tact with  the  wall  throughout  the  instroke.  After  point  "A"  in  Figure  39, 

TcRIT  ^ P ontil  point  "B”  early  in  the  outstroke.  At  this  point,  the  water 
at  this  temperature  would  flash  to  steam.  ~ 

In  Figure  40,  wc  have  employed  equation  (83)  to  plot  the  temperature  gradient 
in  the  water.  The  value  of  B corresponding  to  point  "B”  in  Figure  39  is 
B = 2.35  (t  = .03.3  secs).  Since  2/aO  = .0004882,  the  temperature  rise  above 
the  bulk  (central)  value  is  negligible  at  a distance  of  .00073  ft  (.0088  in) 
from  the  wall.  This  is  comparable  with  the  thickness  of  the  laminar  sublayer 
in  the  boundary  layer,  thus  justifying  our  assumption  of  laminar  flow  values 
for  the  water's  thermal  constants. 


As  an  independent  check  on  tliis  conclusion,  we  note  that  the  example  pulsejct 
consumes  1.65  Btu/pulse.  If  all  this  heat  were  used  to  make  steam,  the 
corresponding  compressed  liquid  volume  would  be  .0383  in^.  The  internal 
diameter  of  tlie  boiler  is  0.6875  in.,  so  the  surface  area  in  contact  with 
the  water  is  10.37  in^.  Spreading  .0383  in^  of  water  over  this  gives  a 
thickness  of  .0057  in,  or  about  half  the  value  computed  (for  the  quasi- 
triangiilar  temperature  distribution)  at  the  interface.  Tlius,  the  two  results 
are  in  reasonable  agreement. 


An  Alternative  Model  for  Water  Heating 

In  the  previous  section,  wc  assumed  that  the  heated  water  moves  at  the  average 
speed  of  the  entire  column  over  the  boiler  surface.  Then  we  found  tliat  its 
thickne.ss  was  less  than  tlic  laminar  sublayer  in  the  boundary  layer,  so  that 
it  is  liardly  moving  at  all!  This  suggests  that  we  should  change  the  model, 
and  assume  tliat  the  advancing  interface  leaves  a thin  film  of  water  "stuck" 
to  the  boiler  wall  ns  it  passes,  so  that  tlic  highest  .surface  water  temperature 
is  (for  this  model)  attained  by  the  water  layer  left  at  the  enti-ance.  And 
so  this  is  also  where  the  wall  temperature  will  be  coolest.  The  concept  of 
"film  resistance"  is  meaningless  with  this  model;  instead,  we  can  compute  the 
water  temperature  gradient  from  the  semi-infinite  slab  result 
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0 0.2  0.4  0.6  O.e  10 

JL  c interface  distance  from  boiler  entrance 

*e  ' MAXIMUM  INTERFACE  PENETRATION 

30.  Water  surface  Temperature  aiul  pressure  at  the  interface  Iq  O^  assuming 
n constant  hoiler  wall  teiiiper.iture  of  .3:’irr'  (h  - .3756  Htu/ft  “hr.  "t') . 
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TIk'  iiont  inli.i'.'d  to  t ho  w.ifor  is  a[>|)roxiiiintcly 


't ii P.p y.Cp f'C^  - icri'o(o)  (90) 

when-  jcrl'oCo)  - 0.96  12  f'l'oiu  oi  ror  function  tables. 

Prosunial)  ly , when  p < ft'f  thin  layer  of  water  i .s  lifteJ  off  tlio  wall  by  a 

layer  of  steam,  ainl  ' starts  to  move  along  the  I^oikr  until  reattaching 

when  p > P(>|,(.,.. 

If  we  put  the  previously  usctl  values  for  the  example  pulsejet  into  equation 
(90)  to  get  a rough  estimate,  we  find 

Q - 0.91  Btu/jnilse,  assuming  t - .0.99  secs 

'I'his  is  about  one  third  of  the  measured  input  of  1.6.9  Btvi/pulse.  I’crhaps  the 
discrepancy  is  due  to  the  fact  that  the  boiler  wall  is  rough,  with  irregulari- 
ties larger  th;m  the  laminar  sublayer,  and  to  heat  conduction  along  the  pipe, 
and  hence  to  the  water  column.  Until  this  is  resolved,  there  is  little  jioint 
in  attempting  to  compute  the  steam  pressuri-t ime  history  on  the  outstroke, 
and  hence  the  boiler  efficiency. 

It  would  also  be  desirable  to  model  the  boiler  and  water  interact  Lon  more 
precisely,  taking  account  of  the  steady  state  lic.at  flow'  to  the  boiler,  and 
the  wall  temp('rature-t  ime  history  as  well  as  that  of  tlu^  water.  The  heat 
balance  integral  tecbnitiues  developed  by  Coodinan^'^  and  Zien^'^  should  probably 
be  tiied  first,  before  facing  the  cost  of  finite  difference  schemes. 


flow  in  a Bleeder  fube 

A blcedci'  tube  is  a small  c.ii)illar>  cfninectcd  tt>  1 be  boiler,  or  to  tlu'  duct 
close  t'l  tlie  boiler,  the  otlu^r  end  bcinj;  free,  or  imniersoJ  in  water.  Its  func- 
tion is  to  bli-ed  gas  from  tlie  boiler,  particularly  when  the  cngiiu  is  operating 
verticallv,  and  so  is  not  sc  I f- purgi ng. 


If  u;'  iss'i::;.'  that  the  flow  in  the  hleecier  tiiho  is  lamiruir,  so  that  the  D'Anry 
frier  i'ln  faetof  f is  jii^’t'ii  h)- 

r.l  64-.) 

^ " Re  " lid 

then  from  c-quat  i on  (7),  tlie  equation  of  nuition  for  the  water  inside  the  tube  is 


I>„  - P 


(:(h  ■■  I 


• • * ' • 

XX  + n(x)  + Kxx  + gx 


where  It  is  the  vt'rtical  distance  hetvvoen  tlte  top  of  the  hleeJcr 
tube  and  its  far  end 

I,  is  the  tube  length. 

Since  the  bleeder  lube  is  long  in  relation  to  its  diameter,  it  seems  re.ason.-'hle 
to  wi'ite  a - (I  ror  inward  motion  as  well  as  outwards.  Also,  the  total  iitoeement 
in  the  tube  is  small  in  relation  to  its  length,  so  we  can  replace  gx  with  gL. 
'I'hei'r  equation  (Ol)  simplifies  to 


P,..  - P 


gh  ^ (x  + K.\)l,  - Pj  say 


= Pj/I.K  + (x^  - Pj/l.K)e'^^ 

and 

X -=  x^^  + Pjt/l.K  - (1/K)fx^^  - lyiKle'  (!>4 ) 

Ihiring  this  iinjnilsive  ]ieriod  of  the  cycle,  the  velocity  is  changed  from  Xj  to 


Xq.  1.0. 


X^  = Xj  - 1/(,1-K  (nro 

where  I is  defined  l>y  equations  (.^4)  and  (35).  Thus  from  equation  foi- 

a cycle  period  of  tp 

Xj  - Pj/i,K  ^ (ij  - 1/pI.K  - Pj/LK)e''^^P 

- Pj/bK  - (l/pl.K)le''^^P/(l  - e''^^P)|  (9(.) 

x^  - Pj/bK  - (l/pl.K)(  1/(1  - (117) 

and  knowing,  the  values  of  the  constants,  one  can  eompute  whether  Xp  - x^  over 
a cycle  is  positive  or  negative. 
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Some  experimental  measurements  of  net  bleeder  flow  are  given  in  Figure  41, 
obtained  with  a vertically  operating  pulsejet.  IVhcn  the  net  flow  was  into 
the  engine,  it  stopped  running  quite  quickly.  When  the  net  flow  was  outwards 
it  ran  steadily  and  pulses  or  bursts  of  bubbles  were  seen  to  issue  from  the 
ends  of  the  bleeders  each  time  the  main  body  of  water  was  reversed  in  the 
boiler. 


» 
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BLEED  TUBE  EXHAUST  HEIGHT  ABOVE  WAl 


FLOW  RATE  IN  OUNCES  PER  STROKE 
(Minus  indicates  a Net  Flow  into  the  pulsejet)  8/9/76 


Figure  41.  Bleeder  Tube  Flow  Rate  as  a Function  of  Tube  Height  Above  Water  Level 
for  Various  Configurations.  Bleeder  tube  I.D.  = O.IS  in. 
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SOMi;  F.XPI-.RIMHNTAL  OBSFRVAT[ONS 


Zll'Ci' ‘tnd  Pressure  Mcasuremon t s 

Figure  42  contains  typical  traces  of  force  and  internal  pressure,  the  latter 
being  measured  at  the  boiler  entrance.  The  pronounced  "ringing"  is  due  to 
the  impulsive  nature  of  the  water  column's  reversal  in  the  boiler. 

The  comljination  of  shock  loading  with  rapid  temperature  change  presents  a 
very  challenging  environment  for  the  experimentalist.  (Ke  have  also  spent  a 
significant  part  of  our  budget  on  the  repair  of  supposedly  rugged  transducers.) 
There  seems  little  hope  of  measuring  average  thrust  from  data  such  as  Figure 
42,  and  we  changed  to  measuring  the  time-dependent  fluid  velocity  near  the 
exit.  One  check  on  the  accuracy  of  such  measurements  was  to  separately 
integrate  the  instroko  and  outstroke  records,  to  see  if  they  agreed  together, 
and  with  the  stroke  actually  observed  in  the  transparent  duct.  Unfortunately, 
this  docs  not  validate  the  thrust  computation,  because  even  a low  frequency 
response  record  will  still  integrate  to  the  correct  stroke  value.  However, 
comparisons  with  our  final  solution  did  show  reasonable  agreement. 

The  final  solution  for  average  thrust  measurement  was  also  the  simplest.  The 
engine  was  mounted  in  a heavy  "boat"  floating  in  a pool,  and  the  average  thrust 
was  taken  to  be  the  "bollard  pull,"  the  latter  being  measured  with  a spring 
balance. 


At  the  higher  boiler  temperatures,  the  beat  required  by  some  of  the  boilers  of 
Figure  9 was  more  than  the  6250  watts  available.  In  such  cases,  the  boiler 
was  heated  to  a high  temperature  before  the  engine  was  started.  Then  the 
engine  was  started,  either  by  inserting  a finger  up  the  tail  pipe  and  with- 
drawing it  rapidly  ("hand  starting")  or  by  squirting  cold  water  in  through 
one  of  the  "bleeders."  Then  thrust  measurements  were  made  as  the  boiler 
temperature  fell,  with  results  like  Figure  43. 

An  indication  of  the  boiler's  thermal  inertia  is  given  by  Figure  44.  IVhen  the 
engine  duct  is  in  the  water,  but  is  not  running,  power  is  switched  off  at 
time  t = 0,  but  the  temperature  at  the  thermocouple  station  continues  to  rise 
for  the  next  minute,  as  the  temperature  gradients  change  to  equilibrium.  Then 
it  starts  to  cool  off,  most  of  the  boat  being  lost  to  the  water,  it  is  presumed, 
by  conduction  along  the  pipe. 


When  the  power  is  not  switched  off,  and  the  engine  is  started  instead,  the 
temperature  falls  rapidly,  as  shown  by  the  lower  curve  in  Figure  44,  until  it 
reaches  a stable  temperature.  The  power  which  would  have  been  required  for 
stable  operation  is  computed  by  differentiating  the  curve,  witli  results  shown 
in  Figure  45.  This  is  based,  of  course,  on  the  relationship 
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AVERAGE  THRUST  IN  POUNDS  CYCLE  FREQUENCY  IN  HERT: 


Date:  H'SOpb 


Confijjiirat  ion : Mill  t i -Hole,  Electrical 
Heaters 

Length:  5'  0.S7S"  T.D. 

Angle:  4®  Boiler  Down, 

Unit  has  (2)  reed  Water  Tubes  Near  Boiler 


Figure  43.  Average  Thrust  and  Cycle  Frequency  as  a Function  of  Boiler  Temperature. 
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where  is  the  boiler  woir.ht  and  (ip  its  specific  heat.  Figure  46  shows  that 
riie  response  tinic>  of  the  thei'i:iocouplc  is  short  enoiigli  for  this  methodology  to 
be  reasonably  meaningful. 

Figure  47  shows  the  effect  of  duct  angle  on  thrust  and  power  rcc|uiremcnt  for 
the  single  hole  boiler  coated  witli  scale.  Continuous  operation  was  possible 
for  duct  angles  between  -3°  (boiler  liigh)  and  +4°  when  a bleeder  was  installed. 
The  lower  limit  is  close  to  the  point  where  stripped  boundary  layer  water 
can  run  into  the  boiler.  The  upper  limit  is  perhap.s  connected  witli  the 
engine's  ability  to  purge  itself  of  gas  so  that  it  can  escape  around  the  duct 
bend  into  the  water  column.  Yet  a month  or  two  earlier,  the  same  engine 
ran  vertically  (a  = -90“ ) for  substantial  periods!  Possibly,  some  leak 
developed  in  the  boiler,  which  only  opened  when  the  boiler  was  hot,  and 
therefore  defied  detection.  Alternatively,  as  will  he  mentioned  later,  boiler 
wall  scaling  may  have  a pronounced  effect  on  engine  performance.  Pcriiaps  the 
engine  operated  in  vertical  mode  while  the  boiler  was  clean.  Certainly,  it 
was  considerably  scaled  in  the  later  tests. 

Figures  48  and  49  show  the  effect  of  some  configuration  changes.  One  of  the 
configurations  - the  "phase  valve"  is  illustrated  in  Figures  50  and  51.  As 
can  be  seen,  the  valve  opens  on  the  insti'oke,  and  it  was  hoped  tliat  this  would 
force  a significant  quantity  of  cold  water  into  the  steam  field  and  promote 
rapid  condensation.  In  practice,  it  had  little  effect  on  the  engine's  opera- 
tion, probably  because  the  water  flow  was  negligible  compared  with  the  volume 
of  the  stripped  boundary  layer  water. 

Experiments  with  a steady  stream  of  cold  water  injected  through  the  nozzle  from 
an  external  faucet  showed  a very  marked  effect,  presumably  because  the  higher 
pres.sure  gave  a much  larger  flow.  As  the  faucet  was  opened,  the  stroke 
decreased  and  the  frequency  increased.  At  high  flow  rates,  the  unit  went  into 
"fibrillation,"  the  column  moving  with  very  small  amplitude  and  very  high 
frequency,  with  no  discernible  thrust.  Then,  as  the  faucet  was  closed,  it 
returned  to  normal  operation. 

It's  hoped  that  if  such  high  pressure  injection  could  be  confined  to  the  in- 
.stroke,  more  powerful  operation  (more  thrust)  could  be  achieved.  We  hope  this 
because,  on  "start-up,"  a pulsejet  delivcr.s  about  twice  as  much  thrust  as  normal, 
and  operates  at  a higher  frequency.  Three  factors  could  account  for  this  - 
colder  stripped  boundary  layer  water,  less  trapped  gas,  and  a hotter  boiler 
wall  surface.  Phased  injection  of  cold  water  sjiould  improve  boiler  penetration 
and,  hopefully,  improve  trapped  gas  purging,  and  it  should  increase  the  fre- 
quency by  reducing  the  duration  of  the  in.stroke. 

During  the  experiments,  a bifurcated  engine  was  operated  with  a head  of  up  to 
six  feet  on  its  inlet  leg,  apparently  without  influencing  its  operation.  The 
frequency  and  stroke  was  unaffected  by  this  head,  so  that  the  thrust  (which 
was  not  measured  directly)  was  presumably  unchanged. 
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Thermocouple  in  Air  80°F 
Thermocouple  in  Water  75°F 
Thermocouple  in  Copper  Boiler  SO 
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Figure  46.  Temperature  Response  of  the  Water  Pulsejet  Thermocouple  in  Various  Media.  (8/4/76) 
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Static  Thrust  for  Various  Angles  of  Duct.  Single  hole  boiler,  5"  long 
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Figure  49.  Boiler  Temperature  vs.  Static  Thrust  for  Various  Pulsciet  Configurations.  Single  h 
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Por  coiv.jiarat  i VC 
familiar  form  of 
hydrocarbon  foci 


purposes,  the  data  presented  can 
"specific  fuel  consumption"  by 
per  hour  (Wp) , i.c. 

F 20, non 


(Ib/hr) 


be  converted  to  a more 
relating  KK  to  pounds  of 


The?)  for  the  multi-hole  boiler  engine  of  Figures  4;i  and  IS 

Boiler  temperature  = 270  300  350  400  500  600  °F 

s.f.c  (-  Kp/T)  = 2.06  1.95  1.45  1.33  1.33  1 . 3S  Ib/lb.hr. 

(As  me)itio))cd  below,  the  lowest  figtire  acliieved  during  this  investigation  was 

0.45  Ib/lb.hr.  with  a cleaned  version  of  the  same  boiler.)  To  obtain  overall 
sfc,  these  figures  must  be  divided  by  combustion  and  boiler  efficiencies. 


Boi 1 er  Seal ing 

A thin  deposit  of  calcium  forms  on  the  boiler  wall^*^  after  a few  hours  of 
running.  Scrapings  of  this  have  the  following  relative  volumes: 


Calcium 

I00"i 

(probably 

carbonate) 

Copper 

10% 

(possibly 

an  oxide) 

Zinc 

3% 

1 ron 

1% 

head,  strontium  and 

manganese 

0.1% 

In  one  series  of  experiments,  the  performance  of  a multi-hole  boiler  engine 
was  measured  before  and  after  mechanically  reaming  the  boiler  to  remove  the 
scale.  As  shown  in  Figures  52  - 59  and  Table  4,  its  performance  after  reaming 
was  significantly  improved.  This  may  indicate  that  the  thermal  resistance  of 
the  scale  is  detrimental.  Alternatively,  the  reaming  may  have  made  the  wall 
surface  substantially  rougher,  which  would  also  have  increased  the  heat  trans- 
fer. It  '.iras  not  possible  to  follow  up  on  this  boiler  to  resolve  the  question, 
because  it  developed  a crack  which  put  an  end  to  its  operational  life. 

For  these  experiments,  the  water  pulsejet  was  equipped  with  a multi-hole  boiler 
and  a bell  mouth  nozzle.  The  angle  of  duct  incline  varied  from  +6°  to  +12®, 
depending  on  duct  length.  The  thrust  measuring  devices  included  a total  head- 
static  system  at  the  end  of  the  pulsejet  duct  and  a thrust  plate  mounted  in 
the  tank,  six  inches  from  the  nozzle. 

The  total  head-static  tubes  were  connected  to  a differential  pressure  trans- 
ducer whose  output  was  recorded  on  a Honeywell  Visicorder.  A thrust  plate 
consisting  of  a 17-1/2"  x 1"  x 3/8"  steel  beam  with  strain  gauges  attached, 
and  a 7-7/8  inch  diameter  aluminum  plate  mounted  on  the  bottom  of  the  beam 
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AVERAGE  THRUST  IN  POUNDS 


CYCLE  TIME  IN  SECONDS 


AVERAGE  THRUST  IN  POUNDS 
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/,  CYCLE  TIME  G THRUST ! AFTER  BOILER  REAMED  OUT 
■"CYCLE  TIME  ) THRUST ! BEFORE  BOILER  REAMED  OUT] 
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Figure  !>3.  Runs  290-298,  duct  length  = 7.2.')  foot. 
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CYCLE  TIME  IN  SECONDS 


AVERAGE  THRUST  IN  POUNDS 


CYCLE  TIME  IN  SECONDS 


AVERAGE  THRUST  IN  POUNDS 


3.0 


/ CYCLE  TIME  OTHRUST!  AFTER  BOILER  REAMED  OUT 
/''CYCLE  TIME  □ THRUST ! BEFORE  BOILER  REAMED  OUT 


CYCLE  TIME  IN  SECONDS 


AVERAGE  THRUST  IN  POUNDS 


Figure  S6.  Runs  310-315,  349-354,  duct  length  = 5.75  feet. 
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CYCLE  TIME  IN  SECONDS 


CYCLE  TIME  IN  SECONDS 


AVERAGE  THRUST  IN  POUNDS 


AVERAGE  THRUST  IN  POUNDS 
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BOILER  TEMPERATURE  IN  "F 


l igiirp  59.  Runs  327-552,  30)6-369,  iliict  length  ~ 4.25  I'oet . 
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was  list'd  111  sa;;n.'  of  the  tests.  It  was  calihratoJ  by  c lamp i ag  the  steel  be.am 
in  a cantiK'Va’.'  ftishioii  and  li.'ingiiig  biiowii  ivcights  from  it.  The  output  f j'oiii 
tiu’  sfiain  g.inyos  was  .also  tecorded  on  the  V'i  si  condor . fixcopt  for  tlie  runs 
m.irkcd  "start-up",  dat.a  was  t.akcn  only  after  all  temperat  lues  had  .stabilised. 

The  thrust  figures  recorded  with  the  velocity  head  before  re.iming  sliowcd 
generally  g(>od  agreement  with  the  "flotiting  boat"  data,  which  we  take  to  be 
umpiest  i onribly  correct  because  of  its  simplicity. 

As  sliown  by  figures  60  and  61,  there  is  some  trend  to  lower  thrust,  at  a given 
power  level  , as  the  length  of  the  duct  is  decreased.  And,  in  the  shorter  duct 
Icng.ths.  and  lovv'er  jiowcr  settings,  cleaning  the  hoiler  docs  not  have  much 
effect  on  performance.  The  two  longest  engines  were  improved  markedly  by 
boiler  cleaning,  csiiccially  at  th('  highest  temperatures,  where  their  thrust 
was  doubled.  The  best  "steady  state"  thrust  of  2.. '^5  Ih  at  62.S0  kiV  corresponds, 
on  the  criteria  used  earlier,  to  a si)ecific  fuel  consumption  of  only  O.d.'i 
Ih/lb. hr. 

Figure  62  presents  smootlied  curves  (cross-plot  ted  twice)  of  the  pulsejct 
operating  frequency.  With  the  cle.an  boiler,  it  can  be  seen  that  both  duct 
length  and  boiler  temperature  influence  the  cycle  time;  a.s  is  generally  the 
case,  whatever  the  condition  of  the  boiler.  (This  boiler  has  also  been 
operated  at  -S  llz,  using  a 2.75  ft  duct.) 

For  some  unknown  reason,  the  engine  could  not  be  operated  when  the  boiler 
temperature  was  in  tlic  range  375-475°F. 


Thrust  Plate  Experiments 

During  and  immediately  after  World  War  II,  the  thrust  of  air-breathing  pulscjets 
was  often  measured  by  a "thrust  plate"  placed  in  their  e.xhaust.  This  appears 
to  be  unfeasible  in  water,  as  Figure  63  attests.  "Ringing"  of  the  plate  is 
quite  severe,  but  most  of  the  error  is  thought  to  be  due  to  the  large  virtual 
water  mass  associated  with  the  plate,  and  the  large  scale  vorticity  generated 
in  the  tank  by  the  successive  pulses  from  the  engine. 
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HEAT  INPUT  IN  WATTS 
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liuilcf  u i t li  si'.'ili'.  riic  iiuiiiIk’i-s  I'.ivu  thu  durt  liuiitth 


-lioh' 

1 1)  fl'Ot  . 
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HEAT  INPUT  IN  WATTS 
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CYCLE  TIME  IN  SECONDS 


THRUST  FROM  THRUST  PLATE 


CONCl.llSIONS 


1.  TIu-  McHugh  cycle  is  qurint  i tat  t ve  ly  cxplavnahle  by  the  well-known  laws  of 
thormoJ)  tiaini  cs , although  much  remains  to  bo  done  before  an  adequate  pre- 
dictive ca])ability  for  design  purposes  can  be  claimed. 

2.  Kliile  it  docs  not  yet  approach  the  efficiency  of  a conventional  steam 
engine,  cycle  efficiency  may  be  adequate  for  some  applications  where 

cost  is  more  important  than  fuel  consumption.  The  best  measured  efficiency 
corresponds  to  a specific  fuel  consumption  of  Ib/lb.hr  (assuming  80^ 

b<nlcr  efficiency)  which  is  about  five  times  the  fuel  consumption  c>f  a 
diesel  engine  driving  a water  propeller  producing  8 Ib/hhp.  Since  pulst>- 
jet  technology  is  i7i  its  infancy,  this  gap  cat)  presum.ihly  be  narrowed. 

7).  The  imlsejet  has  been  shown  to  convey  heat  along  an  0.875  inch  I.D.  duct, 
up  to  7.25  feet  in  length,  at  a rate  of  10  KIV.  There  i s no  indication 
that  this  is  an  upper  limit  for  either  lieat  flow  rate  or  for  length. 

This  performance  compares  well  with  conventional  water  filled  heat  pipes, 
which  typically  choke  at  much  lower  power  levels. 
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